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Zur quantitativen Bestimmung kleiner Zusätze 
zu Alkalihalogenidkristallen 


Von E. Mollwo und R.W. Pohl 
(Mit 1 Abbildung) 


E. Madelungs Vorliebe für Probleme des Kristallbaues sowie 
seine Erfolge auf diesem Gebiete sind bekannt. Trägt doch z. B. eine, 
für das Verständnis des Ionenkristallbaues wichtige Konstante seinen 
Namen. Es erscheint uns daher nicht unangebracht, in diesem Fest- _ 
heft über eine Arbeit zu berichten, die Änderungen im Kristallbau 
eines Alkalihalogenids behandelt. Der Zweck dieser Untersuchung ist 
es, ein neues Eichverfahren für die optische Mengenbestimmung kleiner 
Zusätze zum Kristallmaterial zu entwickeln. 

Seit vielen Jahren werden in unserem Institut photochemische 
Umsetzungen sowie Fragen der Elektronenleitung in ee 
kristallen untersucht. Stets spielen dabei gewisse äußerst geringe 
Zusätze zum reinen Kristallmaterial eine entscheidende Rolle. Für 
ihren Nachweis ist wegen der kleinen Konzentration der einfachste, 
oft allein mögliche Weg die Ausmessung des optischen Pre 
spektrums. Handelt es sich dabei um eine einfache glockenförmige 
Absorptionsbande, so kann man nach einer Formel der Dispersions- 
theorie auch eine quantitative Bestimmung vornehmen. Die Formel 
lautet 


1) H. 


Dabei bedeutet N die Zahl der Absorptionszentren pro Kubikzenti- 
meter, A „ax die Absorptionskonstante im Bandenmaximum gemessen in 
em~?, H die Halbwertsbreite in e-Volt, und A eine aus der Dispersions- 
theorie folgende Konstante der Größenordnung 1016 em-*-e-Volt-!, 

Leider gehen in die Herleitung dieser Konstanten 4 lies 
kaum nachprüfbare Annahmen über den Absorptionsmechanismus ein. 
Daher bleibt eine Mengenbestimmung auf diesem Wege stets mehr 
oder weniger unsicher. Um die einfache optische Methode auch 
quantitativ ausnutzen zu können, bleibt daher nur eine Fichung 
mit irgendeinem anderen Verfahren übrig. 

Als solche Eichverfahren sind bisher von uns zwei ben Beide 


max 


angewendet worden, ein elektrisches und ein chemisches. Beide 
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wurden benutzt zur Messung der Zahl der Farbzentren, also in Alkali. | 
halogenidkristallen überschüssiger Alkaliatome. Das elektrische Ver. 
fahren beruhte auf der Ausmessung der Stromzeitkurve bein: Aus- 
wandern der Farbzentren aus dem Kristall im elektrischen Felde. 
Das Ergebnis der Messungen war recht unsicher, da der Leitungs- 
mechanismus noch nicht genügend geklärt wart), Man fand für das 
Verhältnis der elektrisch und optisch gemessenen Mengen 
N ektrisch 
N optisch 
Beim chemischen Verfahren wurden die farbzentrenhaltigen Kristalle 
in Wasser aufgelöst, die Alkalität der Lösung nach einem empfind- 
lichen kolorimetrischen Verfahren bestimmt und daraus die Zahl der 
Alkaliatome berechnet. Dabei ergab sich beispielsweise für den Fall 


Werte zwischen 1 und 


der Farbzentren in KC]?) 

Vchemisch — 1.24. 

Noptisch 
Der aus der Dispersionstheorie folgende Wert der Konstanten 4 für 
die Farbzentrenbande in KC] beträgt 1,06 em~-?-e-Volt-!. 

Als drittes haben wir jetzt ein Eichverfahren ausgearbeitet, das 

man als gasvolumetrisches bezeichnen kann. Es beruht auf folgender 
Beobachtung?): Zahlreiche Zusätze zum Schmelzfluß von Alkali- 
halogenidkristallen lassen sich in geringer Konzentration in den festen 
Kristall einbauen. Läßt man in einen solchen Kristall bei hoher Tem- 
peratur einen weiteren Stoff hineindiffundieren, so beobachtet man 
oft chemische Umsetzungen im festen Körper?). Dabei entstehen Im 
unter Umständen neue Stoffe, die sich nicht in dem Maße im Kristall- 
material auflösen, wie es die Ausgangssubstanzen taten. Dann tritt 
eine Ausscheidung ein. Ist der auszuscheidende Stoff ein Gas, s0 tb. 
bildet sich im Kristall eine große Anzahl kleiner kubischer Hohlräume. 


Ab 


Sie liegen mit ihren Flächen parallel den Spaltflächen des Kristall. 2 
Das ganze Kristallstück vergrößert dadurch sein Volumen oft um ‘iin 
mehrere Prozent. In den kleinen Hohlkristallen befindet sich das aus- steh 
geschiedene Gas unter dem Druck einiger Atmosphären. Das sieht man Gen 


leicht beim Auflösen des Kristalls in Wasser. Es entsteht dabei aus 
jedem Hohlkristall eine Gasblase. Das Volumenverhältnis läßt sieh 
bequem mit dem Mikroskop ausmessen und liefert direkt den Druck 
im Inneren der Hohlkristalle, Abb. 1. 


1) E. Mollwo u. W. Roos, Gött. Nachr. Math.-Phys. Kl. (N. F.) 1. S. 10. al 


Die 


1934. 
2) F.G. Kleinschrod, Ann. d. Phys. [5] 27. S. 97. 1936. 
3) E. Mollwo, Gött. Nachr. Math.-Phys. Kl. 8. 51. 1941. « 
4) z.B. Sait Akpinar, Ann. d. Phys. [5] 37. ge 1940. 
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Un: die auf den Gesamtkristall bezogene Konzentration des aus- 
geschiedienen Gases zu bestimmen, braucht man also nur das Ver- 
hältnis Gesamtvolumen der Hohlkristalle zum Gesamt kristallvolumen 
zu messen, und mit dem Druck im Innern der Hohlkristalle zu redu- 
zieren. Wegen der geringen Größe der Hohlkristalle, Kantenlänge 


a 


4 ®_ 


D 
D 

* 
Abb. 1. Beim Auflösen des Kristallstückes entsteht aus jedem Hohlkristall cine 


Gasblase. Lösungsmittel: Glyzerin — Wasser, 10: ] 


etwa 10-2 mm, lassen sich so noch sehr kleine Gaskonzentrationen 
genau bestimmen. 

Mit Hilfe dieses Verfahrens hat jetzt Knoop!) eine Eichung der 
optischen Konzentrationsbestimmung ausgeführt, und zwar für einen 
Zusatz von K,CO, zu KCI. Durch diesen Zusatz erhält der KUl-Kristall 
im ultravioletten Spektralgebiet eine glockenförmige Absorptions- 
bande mit einem Maximum bei 204 ma. Wir haben sie nach früheren | 
Absorptionsmessungen?), sowie nach photochemischen Unter 
suchungen?) der K-O-Bindung zuzuschreiben. Läßt man in 
solehen Kristall bei hoher Temperatur Chlor hineindiffundieren, 
verschwindet die Absorptionsbande praktisch vollständig. und es ent- 
stehen die oben erwähnten Hohlkristalle. Irgendeine neue Absorp- 
tionsbande war in dem Spektralbereich 2 > 190 my nicht festzustellen. 
Die Umsetzung lautet in chemischer Schreibweise: 


2K,CO, + 2Cl, = 4KCl + 200, + 0,. 


s ausgeschiedene Gas ist also ein Gemisch von CO, und O,, und 
r entsprechen drei ausgeschiedenen Gasmolekülen (2CO, + O,) 


1) E. Knoop, erscheint demnächst. 
2) A. D. v. Lüpke, Ann. d. Phys. [5] 21. S. 1. 1934. 
3) E. Mollwo, Gött. Nachr. Math.-Phys. Kl. (N. F.) 3. S. 149. 1939. 
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zwei im Kristall eingebaute Zusatzmoleküle (2K,CO,). Man muß 
2/ 


also die gemessene Gasmolekülkonzentration mit ?/; multiplizieren, um 
die Konzentration der umgewandelten Zusatzmole küle zu erhalten. 
In zahlreichen Versuchen, wegen der Einzelheiten verweisen wir 

auf die Dissertation von Knoop, wurde folgendes festgestellt: Die mit 
Hilfe der Hohlkristalle bestimmte Gaskonzentration war naclı voll- 
ständiger Umsetzung bei gleichem Kristallmaterial stets dieselbe. Es 
kam nicht darauf an, in welchem Chlordruck oder wie lange erhitzt 
wurde. Demnach ist die Löslichkeit der ausgeschiedenen Gase außer- 
ordentlich klein. Sie kann vernachlässigt werden. Das war auch nach 
früheren Versuchen zu erwarten!). Quantitativ ergaben sich für ein 
 ausgeschiedenes Gasmolekül oder 0,67 im Kristall eingebaute K,C0,- 
5 Moleküle 1,02 optisch gemessene Absorptionszentren [Formel (1)]. 
Dieser Wert wurde in dem Konzentrationsbereich von 3,5-10!® bis 
71017 K,CO,-Molekiilen pro Kubikzentimeter mit Abweichungen 


Tabelle 


Mittleres 
Volumen 
in em? 


Anzahl 
pro cm? Druck 
in den Gas- K,CO,- |K-O-Bin- 
Hohl- | molekiile Moleküle dungen 
kristallen pro em? | pro em? | pro em? 
in Atm. 


pro cm? 


der Hohlkristalle 


Hohlkristallvol. 
Gesamtkristallvol. 


= 
= 
D 
= 
= 
= 
= 
RQ 
= 
~ 
~ 


Absorptionszentren 


8,23. 1016 4,18. 106 2,65-10=10 1,11-107=3 3,0 8,32-1016 5,55-1016 111-10" 


meist unter 5°/, gemessen. Als Beispiel geben wir in Tabelle 1 eine voll- 
‘ständige Meßreihe. Die für die optische Bestimmung nötige Kon- 
stante A hat für die Absorptionsbande des K,CO, in KCI die Größe 
0,94 cm~?-e-Volt-?. 
Man findet also nach dem optischen Verfahren mehr Absorptions- 
_ zentren als K,CO,-Moleküle nach dem gasvolumetrischen Verfahren. 
Diese zunächst überraschende Tatsache erfährt ihre einfache Er- 
klärung durch die zwei K-O-Bindungen, die man in der Struktur- 
formel der Chemiker dem K,CO,-Molekül zuschreibt. 


Dann entfallen auf ein ausgeschiedenes Gasmolekül 2 x ?/, = 1,33 
K-O-Bindungen im Kristall. Diese erscheinen optisch als 1,02 Ab- 
sorptionszentren. Die Zahl für das Verhältnis 

Nzasvolumetrisch = 131 
N optisch 


1) W. Honrath, Ann. d. Phys. [5] 29. S. 421. 1937. 
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steht in guter Übereinstimmung mit dem in der Dissertation von 
Kleinschrod für Farbzentren bestimmten Wert. 

Den Grund für den Unterschied zwischen den chemischen Mengen- 
bestimmungen und den optischen möchten wir im wesentlichen darauf 
zurückführen, daß bei der Herleitung der Konstanten A in Formel (1) 
der Einfluß der Nachbarladungen auf das Absorptionszentrum durch 
den Entelektrisierungsfaktor !/, zu summarisch berücksichtigt wird. 
Eine genauere Begründung scheint uns indessen noch verfrüht zu sein, 


Zusammenfassung 

Das von uns seit Jahren angewendete Verfahren, winzige Zusätze 
zu Alkalihalogenidkristallen aus Höhe und Breite der Absorptions- 
bande zu bestimmen, verlangt die empirische Bestimmung einer Kon- 
stanten. Diese wird auf einem dritten Wege in guter Übereinstimmung 
mit den früheren Eichverfahren erhalten. Das neue Verfahren benutzt 
die Ausmessung kleiner Hohlkristalle, die durch Ausscheidung gas- 
förmiger Reaktionsprodukte im Kristall entstehen. 


Göttingen, I. Physikalisches Institut der Universität, April 1941. 


(Eingegangen 12. April 1941) 
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Beeinflussung von Diffusionen in Membranen 
Von Raphael Ed. Liesegang 


(Institut für Kolloidforschung und Kaiser-Wilhelm-Institut für Biophysik, 
Frankfurt/Main) 


Versuche, die sich als Ziel die Aufklärung einiger Austausch- 
vorgänge im Lebenden gesetzt hatten, ließen bald erkennen, daß 
zunächst einmal eine Beschränkung auf möglichst einfache Verhält- 
nisse notwendig sei. Denn die Schnelligkeit des Durchdringens von 
 Methylenblaulösung durch ein darauf schwimmendes Papier (Post- 
kartenstärke und dicker, gut geleimt) zeigte so überraschend große 
Streuungen, daß ein sehr wesentlicher Eintluß äußerer Verhältnisse 
angenommen werden mußte. Ein halbjähriges Suchen nach den Ur- 

sachen der Streuungen brachte zugleich Aufschlüsse darüber, in 
_ welcher Weise unter diesen Verhältnissen die Flüssigkeitsbewegungen 
in solchen einfachen Membranen stattfinden. 


Obgleich für einen Theoretiker niedergeschrieben, blieb die 
Erzählung zahlenfrei. Madelung hat immer wieder gezeigt, daß ihn 
auch technische Aufgaben anlocken. 


Der Schlüssellocheffekt veranlaßte anfangs viele Irreführungen, 
unter denen die Andeutung einer verschiedenen Durchlässigkeit des 
Papiers je nachdem die Unterseite oder die Oberseite aufgelegt war, 
als Gegenstück zum Reideffekt (= einseitige Durchlässigkeit der 
lebenden Froschhaut) besonders zur Verfolgung der Vorgänge an- 
gelockt hatte. 


Die 18 x 24cm großen Schalen mit Methylenblaulésung und dem 
darauf schwimmenden Papier (125—250 g je Quadratmeter schwer) 
standen zur ständigen Beobachtung dicht an einem Fenster. Häufig 
_kam es vor, daß die dem Fenster zugewandte Oberseite des Papiers 
noch farblos war, während 20cm weiter auf der abgewandten Seite 
sich schon ein tiefes Blau zeigte. Am stärksten zeigte sich der 
Färbungsunterschied, wenn der Wind auf dem Fenster stand. Der 
_ Verdacht, daß ein leichter Luftzug das Durchdringen des Farbstofis 
durch das Papier in so hohem Maße hemmen könne, bestätigte sich, 
wenn die Schale mit einer Ecke in einem Abstand von 10 cm von 
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einem Schlüsselloch angebracht wurde. Diese Ecke blieb beim Papier 
noch farblos, während in größerem Abstand die Färbung immer mehr 
zunahm. Noch weit größer wurden die Färbungsunterschiede, wenn 
über das Papier ein beiderseitig lackiertes Papier mit Ausschnitten 
gelegt wurde: Tiefes Blau an den bedeckten Stellen, Farblosigkeit 
an den unbedeckten (wenigstens in Fensternähe). 

Bei dünnen (80—100g) Papieren äußert sich die Zugwirkung 
viel schwächer, weil sie wegen zu raschen Durchschlags nicht lang 
genug anhält. Man müßte dazu den Zug künstlich verstärken. 

Kommt es darauf an, den Durchschlag einer Farbstofflösung 
durch ein Papier zahlenmäßig zu verfolgen, so muB jede Luft- 
bewegung sorgsam ausgeschaltet werden. 

Eine Randwirkung. Läßt man bei den Versuchen mit der auf- 
gelegten Schablone aus lackiertem Papier den Farbstoff so lange 
einziehen, daß sich auch die unbedeckten Stellen etwas anfärben, 
so zeichnen sich die Ränder der Schablone durch farblose Säume 
aus. Das durch die Schnitte freigelegte Innere des Papiers saugt 
Feuchtigkeit an und wirkt so wie Zug. Bedeckt man das Ganze 
nochmals mit einem zusammenhängenden lackierten Blatt. so daß 
der Luftzug ausgeschaltet ist, so bilden sich nur die hellen Randzonen 
aus. War die Schablone kurz vorher mehrfach benutzt worden, so 
bleibt die Randwirkung aus, weil das Papier innen Feuchtigkeit 
aufgenommen hat. 

Knickwirkung. Es war erwartet worden, daß ein vorher ange- 
brachter Knick in dickerem Papier ein kapillares, deshalb sehr 
rasches Durchdringen des Farbstoffs bringen würde. Das Um- 
gekehrte war der Fall: Diese Knicke erhielten sich noch farblos, wenn 
die Umgebung schon tief blau war. Die erste Annahme, daß die 
Luft im Knick die Diffusion verhindere, traf nicht zu. Ursache ist 
vielmehr die durch den Knick erleichterte Verdunstung. 


Eine Umkehrerscheinung kann auftreten, wenn man beim Ver- 
such mit der aufgelegten Schablone den Farbstoff viel länger ein- 


wer) ziehen läßt. Zuerst kam die Vermutung, der Farbstoff’ kristallisiere 
iufig an den freien Stellen aus (Gegenstiick zur Bildung der Salpeter- 
piers lagerstätten) oder polymerisiere dort. Beim Schwimmenlassen auf 
Seite einer Lösung von Ammoniumbichromat (das viel rascher durchschlägt 
der als Methylenblau) trat aber oben kein Auskristallisieren ein. Die 
Der Deutung ist vielmehr: An den freiliegenden Stellen konzentriert 
stofis sich die Farbstofflösung innerhalb des Papiers durch Verdunstung. 


sich, Dringt sie endlich nach oben, so muß sie eine tiefere Färbung 
von veranlassen. 
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Gedanken an osmotische Wirkungen waren durch Vermuiungen 
über den Stoffaustausch im Lebenden nahegelegt worden. Sie 
schienen sich zuerst zu bestätigen. Wird das Papier vorlıer be. 
strichen mit Lösungen von KNO, oder anderen Salzen oder von 
Nichtelektrolyten wie Harnstoff, dann getrocknet und dann auf der 
Farbstofflösung schwimmen gelassen, so erfolgt der Durchschlag an 
den bestrichenen Stellen viel schneller. Der Farbstoff wandert also 
rascher in den Orten mit höherem osmotischen Druck vor (wie bei 
der sogenannten „vitalen Resorption“. Aber ein Vorstrich mit 
Wasser ohne vorherige Trocknung wirkt genau so. Die Vermutung, 
das osmotisch Wirksame reiße mit dem Wasser auch den Farbstoff 
hinein, kommt sehr viel weniger in Betracht als die Deutung: Jene 
Aufstriche verzögern die Wasserabgabe durch Verdunstung. 
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Das kann zu Überraschungen führen: Ein Drucksachenumschlag 
mit rotem Briefmarkenstempel wurde mit dieser Seite nach unten 
schwimmen gelassen. Der Druck erschien auf der Oberseite tief 
dunkel. Die Farbe enthält (hygroskopisches) Glyzerin. Das wirkt 
stärker fördernd auf die Diffusion als der Fettgehalt hemmend. 
Fettspaltung im Lebenden macht Glyzerin aus dem sonst wasser- 
abstoBenden Fett frei. 


Wasser und Gelöstes wandern getrennt. Wasser dringt in etwa 
2Min. an die Oberfläche des (125g schweren) Papiers. Das Durch- 
dringen des Methylenblaus (mit seinem viel größeren Molekül — 
bei Malachitgrün geht es rascher) dauert 5 Std. oder mehr. 
Methylenblau diffundiert um so rascher, je wasserreicher die Membran 
ist. Deshalb die Verlangsamung beim Verdunsten des Wassers. 


Wird eine stark quellfähige Membran hergestellt durch Tränkung 
von Filtrierpapier mit Gelatine, so wird das Durchdringen des Wassers 
durch den Sog der Quellung beschleunigt. Auch der Farbstoff wird 
von der (Juellung mitgerissen. Gerbung der Gelatineschicht mit 
Chromalaun oder Formaldehyd vermindert Quellung und Farbdurch- 
schlag. Noch viel stärker wirkt Tränkung mit Bichromat und Be- 
lichtung, (darauf Auswaschen) auch dann, wenn die anderen Ger- 
bungen vorhergegangen sind. Diese Wirkung beruht zum Teil auch 
darauf, daß hier die Wirkung des Luftzugs, die Verdunstung, er- 
leichtert ist. Das Verfahren ist also auch geeignet zur Prüfung der 
Gerbfähigkeit. Auch beim Bichromat zeigt sich zuweilen ein Um- 
kehreffekt: War die aufgestrichene Lösung an einzelnen Stellen 
auskristallisiert, so schlägt hier der Farbstoff fast mit der 
Schnelligkeit kapillaren Durchdringens durch. Der Zusammenhalt 
der Gelatineschicht ist mechanisch gelockert. 
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Ist die zum Durchtränken benutzte Gelatinelösung zu dünn 
(unter 2"/,), so schlägt die Farbstofflésung in wenigen Sekunden 
kapillar durch. Nicht nur die Zeit erlaubt Durchschlag durch 
Kapillarität und durch Diffusion zu unterscheiden. Ein Wasser- 
aufstrich auf ungenügend geleimtes Papier verzögert den Farbstofi- 
durchschlag erheblich durch Verstopfung der Poren, während er beim 


geleimten Papier den Diffusionsweg vergrößert. 


Frankfurta.M. 
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Ein Demonstrationsversuch zur Emanationsdif sion 
Von Adolf Krebs 


(Mit 3 Abbildungen) 


Die Diffusion von Emanation in verschiedenen Gasen ist einer- 
seits wegen ihrer Eigengesetzmäßigkeiten, andererseits im Hinblick 
auf die sich dabei abspielenden Elementarprozesse von Interesse. 
Ihre Demonstration ist deshalb sowohl zur Ver- 
anschaulichung radioaktiver als auch gaskinetischer 
Vorgänge allgemeiner Art von Bedeutung. Die 
modernen, alphastrahlempfindlichen, praktisch nicht 
nachleuchtenden Leuchtphosphore lassen hier be- 
sonders überzeugende Effekte erwarten, um so mehr, 
als durch entsprechende Steigerung des Szintilla- 
tionsleuchtens relativ kurze Adaptationszeiten zu 
erwarten sind. Man präpariert zu dem Versuch 
Glasröhren geeigneter Dimensionen (die Diffusions- 


halbwertsstrecken im Gleichgewichtszustand sind 

für AcEm: 0,54 cm, für ThEm: 1,97 em, für RaEm: 

144 cm) oder — noch einfacher schmale Glas- 

streifen, die in die Diffusionsröhren hineingestellt 

werden, mit Zinksulfid 9S4X. Zur Erzielung 

gleichmäßiger Schichten überzieht man das Glas 

aa. mit einer dünnen Celluloselackschicht, auf die das 

Emanations- Zinksulfid aufgestäubt wird. Bei Rohren empfiehlt 

DiffusionsgefaB es sich, sie bis zum Trocknen auf der Drehbank 

laufen zu lassen. Das Schema der Versuchsanord- 

nung zeigt Abb. 1. Die von einem Präparat geeigneter Stärke in das 

Gefäß G geschleuste Emanation wird in den Raum E eingelassen, 

von wo sie einige Zeit später nach Öffnen des Schiebers $ in das 

Rohr D diffundiert. Am Aufleuchten der Röhrenwände bzw. des 

Glasstreifens kann dann der Diffusionsvorgang in allen Einzelheiten 
verfolgt werden. 

Einem größeren Kreis läßt sich der Versuch bei Anwendung 

genügend großer Emanationsmengen mit Hilfe einer empfindlicher 

Photozelle vorführen. Wir wählten als empfindlichste Photozelle eu 
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Photonenzählrohr nach Rajewsky!), das in letzter Zeit öfters zu 

physikalischen Messungen und Demonstrationen benutzt worden ist?). 

Es besitzt in unserem Falle noch den Vorteil, auf die gesetzmäßiger- 

weise mit dem Leuchten verknüpfte y-Strahlung bzw. z-Strahlung 
N) anzusprechen, die mit der Folgeproduktbildung gegeben ist. Der 
benutzte Photonenzähler bestand aus einem Messingrohr von 4 em 
Länge und 2,6 em Durchmesser. Als Photokathode diente ein Mantel 
aus Reinaluminium im Inneren des Messingrohres, als Zähldraht ein 
0,05 mm starker blanker Platindraht. Die Arbeitsspannung des 
Zählers betrug bei Luftfüllung und einem Druck von 50 mm He 
1180 Volt. Das Fenster bestand aus einer Cellophanfolie von 25 u 
Dieke und war 1,2 x 0,7 cm? groß. Je nach der Stärke der bei dem 
Versuch auftretenden y-Strahlung wurde ohne oder mit allseitigem 
Bleischutz des Zählers gearbeitet. Der Bleischutz war 1.5 cm stark 
und die Strahleneintrittsöffnung 0,25 x 2,0 cm? groß. Besondere 
Vorversuche über das Ansprechen dieses Zählers auf das Szintillations- 
leuchten und Messungen des Ausbeutefaktors ergaben, daß ein Po- 
Präparat von 60 e. s. E. auf einer Fläche von 1 cm? in 2 mm Ab- 
stand von der Zinksulfidschicht des dicht vor dem Fenster aufgestellten 
Leuchtschirmes die Stoßzahl des Zählers um 3,5 Stöße pro Minute 


erhöht. 
ich Bei Versuchen mit RaFm läßt man einfach die Emanation aus 
- E in D hineindiffundieren, und kann zur Ermittlung der Gesetzmäßig- 
ind keiten entweder bei festem Beobachtungsort die Zahl der diffun- 
Mr dierenden Atome als Funktion der Zeit messen | N = f(t) für x = const.| 
as- oder bei konstant vorgegebener Konzentration den jeweiligen Ort in 
ell Abhängigkeit von der Zeit bestimmen [r = f(t) für N = const.)?). 
un Nimmt man nach 10-minütiger Exposition in einer 100000 mst- 
en Radonator-Emanationsatmosphäre den Glasstreifen aus dem Diffu- 
das sionsrohr heraus und führt ihn am Zähler vorbei, so erhält man als 
= Verteilung der diffundierten Atome bzw. ihrer Folgeprodukte den in 
ye Abb. 2 wiedergegebenen Verlauf. 
ane Bei ThEm-Versuchen empfiehlt es sich, das Präparat bis zum 
un Gleichgewichtszustand im Diffusionsrohr zu belassen. Ein Radiothor- 
pi präparat von 1,3 mg Radium Äquivalent gibt nach 30-minütiger Expo- 
pr sition in einem 25 em langen Rohr ein weithin sichtbares, gesetzmäßig 
Ana mit der Entfernung von der Emanationsquelle abfallendes Leuchten. 

1) B. Rajewsky, Phys. Ztschr. 32. S. 121. 1931. 
dung 2) M. Schein u. B. Stoll, Helv. Phys. Acta 7. 8. 672. 1934; M. Dja- 
| tschenko u. W. Sselegenew, C. R. (Dokladv) de l’Ac. d. Se. de VURBS 29. 
chen 8. 19. 1940. 
e ein 3) R. Fürth, Phys. Ztschr. 26. S. 719. 1925 ' 
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Abtasten des Leuchtstreifens mit dem Lichtzähler ergibt dic Ver. 
teilung nach Abb. 3. 

Sie zeigt, daß die Zahl der diffundierten Atome e-funktion- mäßig 

mit dem Abstand von der Emanationsquelle abnimmt und in Über- 
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Abb. 2. Diffusionsversuch (Ra Em) Abb. 3. Diffusionsversuch (Th Em) 


einstimmung mit dem Gleichgewichtszustand und dem Diffusions- 
koeffizienten beim Fortschreiten um jeweils 2 em auf die Hälfte ab- 
sinkt. Durch Auffangen des aktiven Niederschlages auf geeignet ge- 
formten Elektroden läßt sich der Versuch noch weiter variieren. 


Frankfurt a. M., K.-W.-I. für Biophysik, 3. 4. 1941. 


(Eingegangen 10. April 1941) 


Anmerkung bei der Korrektur: 
5 Während der Drucklegung wurden die Versuche mit einem 10 mg-Ra- 
Präparat weitergeführt; die Leuchterscheinung ist dann so stark, daß ihr 
-gesetzmaBiger Verlauf und die einzelnen Szintillationen bei diffuser Beleuchtung 
des Raumes schon mit bloßem Auge beobachtet werden können. 
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Zur Thermodynamik der Thermodiffusion 


Von J. Meixner 


1. Vorbemerkungen 


Die Thermodiffusion ist von Enskog!') in einer grundlegenden 
Arbeit nach den Methoden der kinetischen Gastheorie behandelt 
worden. Er beschränkt sich von vornherein auf kugelsymmetrische 
Moleküle. Die Ergebnisse sind recht einfach, nur ihre Herleitung 
und die numerische Berechnung der Koeffizienten wie Wärmeleit- 
fähigkeit, Diffusionskonstante usw. ist recht langwierig. Zahlenwerte 
für die Koeffizienten wurden nur für einige spezielle Wechsel- 
wirkungsgesetze bestimmt. Im Hinblick darauf. daß die Versuche 
zur Thermodiffusion nicht nur mit einatomigen Gasen durchgeführt 
werden und daß die ausgewerteten Wechselwirkungsgesetze manch- 
mal auch nicht annähernd den tatsächlichen Verhältnissen ent- 
sprechen dürften, ist es vielleicht nicht unnötig, den Prozeß der 
Thermodiffusion einmal rein phänomenologisch zu betrachten. Man 
begibt sich zwar damit von vornherein der Möglichkeit, irgendwelche 
Koeffizienten zu berechnen, auf der anderen Seite aber zeigt sich 
ein Weg, der verhältnismäßig schnell auf allgemeine Gesetzmäßig- 
keiten führt, die auch im Falle mehratomiger Gase und bei be- 
liebiger Wechselwirkung gültig sind. Eine Einschränkung bezüg- 
lich der Allgemeinheit muß jedoch gemacht werden. Für sehr schnelle 
Vorgänge, bei denen Relaxationserscheinungen auftreten (verzögerte 
Einstellung des thermischen Gleichgewichts zwischen verschiedenen, 
äußeren und inneren, Freiheitsgraden) sowie für den Fall hoher 
Temperaturen mit der Möglichkeit der Dissoziation beim Stoß von 
Molekülen gelten die folgenden Betrachtungen nicht ohne weiteres; 
eine Verallgemeinerung, die solche Erscheinungen umfaßt, ist jedoch 
möglich. 

Die im folgenden durchgeführten Betrachtungen zur Theorie 
der Thermodiffusion scheiden sich in einen thermodynamischen und 
in einen phänomenologischen Teil. Der erste Teil geht von den 
Erhaltungssätzen für die Masse, den Impuls und die Energie aus 
und führt unter Zugrundelegung allgemeiner thermodynamischer 


1) D. Enskog, Kinetische Theorie der Vorgänge in mäßig verdünnten 
(rasen. Dissertation Upsala 1917. 
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Sätze auf den Begriff der Entropieströmung und der lokale: En- 
tropieerzeugung bzw. Energiedissipation; der zweite Teil behundelt 
dann die Ansätze zur Verknüpfung von Teilchenstrom, Energic- und 
Impulsstrom, mit den Größen, die diese Ströme verursachen, näm- 
lich Konzentrationsunterschiede, Temperaturgefälle und Gesch indig- 
keitsgefille. Ein Vergleich mit den allgemeinen Ergebnissen von 
Enskog gibt ein Beispiel für die Gültigkeit des Reziprozitätssatzes, 
den Onsager?) auf Grund des Satzes von der mikroskopischen Rever- 
sibilität mit Hilfe von Schwankungsbetrachtungen aufgestellt hat. 
Dieser Satz gibt eine bemerkenswerte Einschränkung für die 
phänomenologischen Ansätze bei irreversiblen Prozessen. 

Diese thermodynamisch - phänomenologische Methode, ergänzt 
durch den Reziprozitätssatz von Onsager, ist nicht nur auf das 
Problem der Thermodiffusion, sondern auch auf viele andere irre- 
versible Vorgänge anwendbar. So liefert sie ohne weiteres die 
Thomsonsche Beziehung zwischen den Konstanten des Peltier- 
und des Thomsoneffekts!), die Thomson seinerzeit auf rein thermo- 
dynamischem Wege glaubte herleiten zu können?) Ferner gibt sie 
den Zusammenhang zwischen den Konstanten des Ettingshausen- 
effekts und des 1. Ettingshausen - Nernsteffekts mit dem richtigen 
Vorzeichen, 

Die Begriffe der Entropieströmung und der lokalen Entropie- 
erzeugung haben, so zweckmäßig sie für die Behandlung irreversibler 
Prozesse sind, bis jetzt nur vereinzelt Anwendung gefunden, Enskog’) 
hat sie, von der statistischen Deutung der Entropie ausgehend, in 
die kinetische Gastheorie eingeführt. Borelius*) bezeichnet den 
Quotienten aus Wärmeleitfähigkeit und absoluter Temperatur als 
Entropieleitfiihigkeit, ohne diesen Begriff weiter anzuwenden; er 
weist nur darauf hin, daB damit das Wiedemann-Franzsche Ge- 
setz eine besonders einfache Form annimmt. Der Verf. wurde bei 
seiner Behandlung der thermoelektrischen Effekte im Magnetfeld‘) 
sowohl von der thermodynamischen wie von der kinetischen Seite 
her auf diese Begriffe geführt. Bei der thermodynamischen Behand- 
lung stationärer Vorgänge in Mehrphasensystemen hat Meissner‘) 
die mit irreversiblen Strömungsvorgängen und Wärmeleitung ver- 
bundenen Entropieänderungen ausführlich untersucht. Besonders 


1) L. Onsager, Phys. Rev. 37. S. 405. 1931; 38. S. 2265. 1931. 

2) Vgl. J. Meixner, Ann. d. Phys. [5] 35. S. 701. 1939. 

3) D. Enskog, Ztschr. f. Phys. 54. S. 498. 1929. 

4) G. Borelius, Handb. d. Metallphysik I 1. $.371, 379. Leipzig 1935. 
Meixner, a. a. O. 
6) W. Meissner, Ann. d. Phys. [5] 32. S. 115. 1938. 
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zu erwähnen sind die Arbeiten von Jaumann’) und Lohr®, die 
leider an etwas entlegener Stelle veröffentlicht wurden und daher 


wenig bekannt geworden sind. Sie haben in umfangreichen und JE 
allgemeinen Untersuchungen über geschlossene Systeme physikali- u 


scher und chemischer Differentialgesetze zur Formulierung des En- 
tropieprinzips bei irreversiblen Vorgängen die Begriffe der Entropie- 
strémung und der lokalen Entropieerzeugung herangezogen. Ihre 
Methode steht in engem Zusammenhang mit den hier gegebenen 
Überlegungen, abgesehen von der Anwendung des Onsagerschen 


Reziprozitätssatzes. 


2. Erläuterung der Methode; Begriff der Entropieströmung 


Um den einzuschlagenden Weg klar zu machen, soll zuerst 
der Fall der Wärmeleitung in einem festen anisotropen Körper be- 


handelt werden. Sei T die absolute Temperatur, W der Wärme- 
strom, o die Dichte und c, die spezifische Wärme pro Gramm. Der 
Erhaltungssatz für die Energie spricht aus, daß sich in einem Vo- 
lumenelement eine Erwärmung nur auf Kosten zugeführter Wärme 
einstellen kann; es ist daher 


“dt 


(1 06, 


Die Entropie 9S in einem Kubikzentimeter ändert sich dann 


dT I 


infolge der Temperaturerhöhung und es ist wegen dS = ep (von 
thermischer Ausdehnung werde abgesehen) Ft 
18 1 > 
er = divW. 


Diese Gleichung läßt sich auch in folgender Weise schreiben 


> 
3) div = — (W, grad 7) . 


Wird diese Gleichung iiber das ganze wiirmeleitende System 
integriert, das nach außen hin adiabatisch isoliert sein soll, so ver- 


schwindet das Integral über div (W/T), da es nach dem Gaussschen 


1) G. Jaumann, Sitzber. d. Akad. d. Wissensch., Wien, Math.-Naturw. 
Klasse 120. §. 385, 1911; Denkschr. d. Akad. d. Wissensch., Wien, Math.- 
Naturw. Klasse 95. S. 461. 1918. 

2) E. Lohr, Denkschr. d. Akad. d. Wissensch., Wien, Math.-Naturw. 


Klasse 93. S. 339. 1916; Festschr. d. Deutschen Techn. Hochschule in Brünn, 
8.176. 1924, 
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Satz in ein Obertlichenintegral über die Begrenzung des Systems 


> 
verwandelt werden kann, und dort ist W=0. Das Integra! über 
die rechte Seite von (3) ist daher gleich der Entropiezunahme des 
ganzen Systems in der Zeiteinheit und als solche nach dem zweiten 


Hauptsatz der Thermodynamik positiv. Nun weiß man, dab Vv und 
grad T immer einen stumpfen Winkel einschließen (in isotropen 
Medien sind sie speziell entgegengesetzt gerichtet); d. h. die rechte 
Seite von (3) ist nirgends negativ und deshalb liegt es nahe, sie als 
die im Kubikzentimeter pro Sekunde erzeugte Entropie anzusehen 
(sie sei als lokale Entropieerzeugung bezeichnet). Die Entropie- 
zunahme des ganzen Systems pro Sekunde ist dann gleich dem 
Raumintegral der lokalen Entropieerzeugung. Die Entropiezunahme 
eines Kubikzentimeters pro Sekunde braucht jedoch nach (3) nicht 
gleich der lokalen Entropieerzeugung zu sein und das führt dann 


zwangsläufig dazu, WıT als Entropieströmung zu deuten; d.h. wenn 
die beiden erwähnten Größen nicht einander gleich sind, so ist eben 
Entropie zu- oder weggeströmt. 

Damit ist der zweite Hauptsatz der Thermodynamik, der nur 
etwas über die Entropieänderung eines abgeschlossenen Systems 
aussagt, in der Weise verfeinert, daß in jedem Volumenelement in 
der Zeiteinheit eine bestimmte nichtnegative Entropie erzeugt wird; 
sie hängt von den Zustandsgrößen, d.b. bei dem hier vorliegenden 
irreversiblen Prozeß der Wärmeleitung auch von grad T ab. 

Diese Deutung von Gleichung (3) ist deshalb möglich, weil sie 
den Charakter einer Kontinuitätsgleichung für eine Größe hat, die 
nicht erhalten bleibt, sondern eine räumliche Quelldichte besitzt. 
Eine solche Kontinuitätsgleichung für die Entropie läßt sich bei 
allen irreversiblen Prozessen aufstellen, bei denen es überhaupt einen 
Sinn hat, von Entropie und Temperatur zu sprechen, wahrschein- 
lich auch bei solchen Systemen, denen man zwei verschiedene defi- 
nierte Temperaturen am selben Ort zuzuschreiben hat, wie z. B. im 
Gasentladungsplasma (Klektronentemperatur und Temperatur des 
neutralen Gases) und bei der paramagnetischen Relaxation (Spin- 
und Gittertemperatur). 

Dabei erhebt sich die Frage, wieweit es einen Sinn hat, bei 
irreversiblen Prozessen von der Temperatur an einer bestimmten 
Stelle zu sprechen. Eine Temperatur kann die Thermodynamik 
nur für den Fall des thermodynamischen Gleichgewichts definieren. 
Andererseits hat es erfahrungsgemäß einen guten Sinn, wenn man 
in einem wärmedurchströmten Körper jeder Stelle eine bestimmte 
Temperatur zuordnet. Die kinetische Gastheorie kann diese Tat- 
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sache für Gase begründen; sie zeigt zwar, daß im Falle eines be- 
wegten Gases in höherer Näherung eine Erweiterung des Tem- 
peraturbegriffs der Thermodynamik notwendig ist; aber solange die 
Temperaturänderung über eine freie Weglänge klein gegen die abso- 
lute Temperatur ist, kommt man mit dem thermodynamischen Tem- 
peraturbegriff aus; diese Bedingung ist aber, außer in hochverdünnten 
Gasen, praktisch wohl immer erfüllt. In Flüssigkeiten und in festen 
Körpern liegen die Dinge ähnlich. 

Meissner’) zeigt, daB ein sehr dünnes Plättchen mit ideal 
glatter Oberfläche, das sich in einem expandierenden Gas mit- 
bewegt, einmal in die Expansionsrichtung, einmal senkrecht dazu 
gestellt, verschiedene Temperaturen annimmt. Eine Definition der 
Gastemperatur durch das Temperaturgleichgewicht mit dem Plätt- 
chen verliert daher in Strenge ihren Sinn; praktisch sind die Plätt- 
chentemperaturen in den verschiedenen Stellungen im allgemeinen 
so wenig verschieden, daß sich so doch eine brauchbare Temperatur- 
definition ergibt. 

Für die Entropie eines idealen Gases mit Temperaturgetälle 
ohne innere Reibung erhält Enskog?) in zweiter Näherung einen 
Ausdruck der Form 


S = S,(T,p) + «(T,p)-(grad T)?. 


Auch hier kann das Zusatzglied zur thermodynamischen Entro- 
pie S, vernachlässigt werden, solange die Temperaturänderung auf 
einer freien Weglänge klein gegen die absolute Temperatur ist. 
Auch die innere Reibung ändert nichts daran, wenn sich die Ge- 
schwindigkeit über eine freie Weglänge nicht sehr ändert. 

Man kann dementsprechend bei der Wärmeleitung und auch 
bei vielen anderen irreversiblen Prozessen Temperatur und Entropie 
durch die Werte definieren, die sich für beide Größen ergeben, 
wenn man ein hinreichend kleines Volumenelement (das aber 
noch genügend viele Atome oder Moleküle enthält) plötzlich iso- 
liert und das thermodynamische Gleichgewicht sich in ihm ein- 
stellen läßt. 

Der zweite Schritt zur Behandlung des Problems der Wärme- 
leitung besteht nun darin, daß für den Zusammenhang zwischen 
Wärmestrom und Temperaturgefälle ein geeigneter, sogenannter 
phänomenologischer Ansatz gemacht wird. Bei allen Ansätzen dieser 
Art stellt im allgemeinen die lineare Abhängigkeit zwischen den 


1) W. Meissner, Ann. d. Phys. [5] 36. S. 303. 1939. 
2) D. Enskog, Ztschr. f. Phys. 54. S. 498. 1929. 
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3rößen eine gute und oft in sehr weiten Grenzen brauchbare 
 Näherung dar. Speziell für die Wärmeleitung im anisotropen 
Körper ist der allgemeinste solche Ansatz (W, = Komponenten des 


tel 


Die neun Koeffizienten #;, sind die Komponenten des Wärme- 
leitungstensors. Die Zahl der unabhängigen Koeffizienten reduziert 
sich im allgemeinen noch, wenn die Symmetrieverhältnisse des 
anisotropen Körpers berücksichtigt werden. Eine weitere Reduktion 
verlangt im allgemeinen das Prinzip der mikroskopischen Reversibilität, 
bzw. der von Onsager’) daraus abgeleitete Reziprozitätssatz, der 
eine makroskopische Aussage darstellt. Danach muß der Wärme. 
leitungstensor symmetrisch sein (x,, = %,;); er enthält also keinen 
rotatorischen Anteil. Sehr genaue Messungen von Voigt?) und 
Soret°) bestätigen dies. 

Für die lokale I; ergibt sich mit (4) Bee: 


Ausdruck 
3 


ik Ox; Oa, 


Wärmestroms W) 


t=1 k=1 
Wegen ‘+> 0 folgt daraus die bekannte Tatsache, daß der 
u Wärmeleitungstensor positiv definit sein muß. 
Eine allgemeinere Formulierung des Reziprozitätssatzes wird bei 
der Behandlung der Thermodiffusion gegeben. R 


3. Die Kontinuitätsgleichungen für die Thermodiffusion 


- 


Gegeben sei eine Mischung zweier idealer Gase 1 und 2. Die 
Koordinaten ‘seien mit &,,2,, 2,, die Zeit mit t bezeichnet. Die 
Dichten der Komponenten seien 9,, 9, g/em*, die Gesamtdichte 
0=0, + 0,. Die Schwerpunktsgeschwindigkeit eines Volumen- 


2 


elements sei v cm/sec (bezeichnet man mit v, und v, die mittleren 
Geschwindigkeiten der Moleküle 1 und 2 an einem bestimmten Ort, 
> > > > 

so ist also 9,0, + 0,2, =v). D, und D, seien die Diffusions- 
ströme in g cm”? sec~!; für sie gilt 
=> > HE id 


1) L. Onsager, Phys. Rev. 37. S. 405. 1931; 38. S. 2265. 1931. 
2) W. Voigt, Göttinger Nachr. Math.-Phys. Klasse S. 87. 1903. 
) Ch. Soret, Arch. de Geneve 29. 1893. 32. 1894. 
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Die Diffusionsströme geben also an, welche Mengen der beiden 
Komponenten in 1 sec durch eine Fläche von 1 cm? diffundieren, 


die sich mit der Geschwindigkeit v bewegt und auf der Richtung 


von 0 senkrecht steht. 

Auf die beiden Molekülarten mögen pro Gramm die äußeren Kräfte 
K und K, wirken. Die p,, seien die Komponenten des Druck- 
tensors p (Reibungsdruck + Gasdruck p). Wegen des Schwerpunkts- 
und des Drehimpulssatzes für ein abgeschlossenes System (Impuls- 
und Drehimpuls können sich nur durch Einwirkung äußerer Kräfte 
ändern) gilt p,, = p,, fir i,k =1,2,3. Die gesamte Energie von 
1gGas setzt sich aus kinetischer Energie v?/2 und innerer Energie U 
zusammen. Dann gelten die Kontinuitätsgleichungen ') 


: 
+ o,divv + div D, = 0: 


dv 


K, +0, K, — Divp, 


div (p; _(K, D,) 


K,d)-Zu# 


dt 


Die letzte Gleichung ist der Ausdruck dafür, daß die innere 
und die kinetische Energie zunehmen auf Grund der von den 
Kräften an der betrachteten Masse von 1 g geleisteten Arbeit und 


— 
auf Grund der infolge der Wärmeströmung W zugeführten Energie. 
Die Glieder mit der äußeren Kraft kommen so zustande: 8: (K,, v) 
ist die Arbeit, die an o,/e g Molekülen der ersten Art in 1 Sek. 


. . Sint . . 
geleistet wird; wegen o, (v, — v )= D, und wegen (6) ergibt sich 
der in (9) angegebene Ausdruck. 

Durch Addition der gg Gl. w ) entsteht 


10) o div T = (0, 


Ferner sei noch die aus (7) und (10) — Formel 


FAR divD, = 


1) Der Vektor Div p ist die Tensordivergenz von p; v,; = AEDRERRIRE 


> > 
des Vektors r: der Vektor (p: v) hat die Komponenten 


Ss i : 
T=1,2,3. 
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angemerkt. Sie sagt aus, daß sich die Konzentrationen nur infolge 
Diffusion, aber nicht infolge Konvektion ändern. 

Die totale Differentiation bedeutet wie immer die Änderung 
bezüglich der mitbewegten Masse von 1g. Es ist 
d 0 
dt Ot 
Die Kontinuitätsgleichung (9) für die Energie läßt sich noch 


(11) (7, grad). 


umformen, indem (8) skalar mit v multipliziert und der entstehende 


Ausdruck für „ z iB (9) eingesetzt wird. Es ergibt sich 
- at 


D,) + (X, D, )- divW. 


Für äußere Kräfte, die zur Masse proportional sind (Schwer- 
kraft, Zentrifugalkraft), ist K, = K, und wegen (6) verschwindet 
dann der Einfluß der äußeren Kräfte in (12). Eine solche 
Spezialisierung soll jedoch im folgenden nicht angenommen werden, 

Für Mischungen zweier idealer Gase kann man nach der Thermo- 
dynamik die innere Energie und die Entropie additiv aus den ent- 
sprechenden Bestandteilen für die reinen Komponenten zusammen- 


setzen. Es ist 


(13) 


wo u, und u, die inneren Energien der reinen Gase 1 und 2 pro 
Gramm bei der betreffenden Temperatur sind, sie hängen nicht vom 
Druck ab. Die Entropie für 1 g ist 


WO 5; i 8, die Entropien der reinen Gase 1 und 2 pro Gramm 
bei der betreffenden Temperatur und bei den Partialdrucken p, 
und p, sind. Diese errechnen sich aus dem Gesamtdruck p zu 


Für gilt nach der Zustandsgleichung des 
idealen Gases 


(16) p= (+ RT. 


Die M, sind die der beiden Komponenten 
der Mischung, R ist die Gaskonstante. 
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Die Kontinuitätsgleichung für die Entropie läßt sich nun all- 
gemein so gewinnen, daß man (14) total nach der Zeit differentiiert ji 
und die auftretenden totalen Differentialquotienten mittels der 

Kontinuitätsgleichungen (7) bis (9) ersetzt. Es ist y 
S ( a 

Die Entropie eines gewissen Gebiets des Gases, das gerade 1 g 
enthält und mit dem Gas mitströmt, ändert sich also wegen der | 
Änderung der Gaszusammensetzung, der Temperatur und der Partial- 
drucke. Wegen der thermodynamischen Beziehungen 


18) Tds,(T,p,) = du,(T +24(--), 


und wegen (15) und (16) wird 


“2 at | | 
R dg, 


Me d Mo “at” 


Die Glieder, die du, /dt und du,/dt enthalten, können mittels 
der nach ¢ differentiierten Gl. (13) eliminiert werden. Für dU/dt 
t dann der Ausdruck (12) einzusetzen; die Differentialquotienten 
von 0,, 0,, g nach der Zeit sind durch die Gl. (7) und (10) auszu- 
drücken. Dann ergibt sich unter Berücksichtigung von (16) 


= = — 8, + ar) div D, 58. + 


- div W + Pix) 


i=1 k=l 


amm 


) + (K,, D,)]- 


ö,, ist der 
Die Koeffizienten 


d die Planckschen charakteristischen Potentiale Kompo- 

nten der Gasmischung bei der Temperatur T und beim Partial- 
druck p,. Eine kleine Umformung führt nun von (20) auf eine 
Kontinuitategleichung vom Typ (3) 
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mit der Entropieströmung 


G > 
(23) 


und der lokalen Entropieerzeugung 


> 


K, 


+ grad ®, | 


Der Entropiestrom S’ ist demnach nicht allein durch den 
Wirmestrom gegeben; es ist von ihm ein Anteil in Abzug zu bringen, 
der den Gibbsschen thermodynamischen Potentialen (=— T® 
der beiden Gase, multipliziert mit den in 1 Sek. durch 1 cm? diffun- 
dierten Mengen gleich ist. Dasselbe ergibt sich in der Theorie der 


thermoelektrischen Effekte’). In der Entropieströmung S’ ist d 
durch Konvektion bedingte Anteil ebensowenig enthalten, wie 


> 
dem oben eingeführten Wärmestrom W. Er ergibt sich aus (22), 
wenn man die Entropieänderung eines im Raume festen Volumens 
von 1cm* berechnet. Diese ist nach (10), (11) und (22) 


6 (0 S) 
Ot 


= div (vo S48"). 


(25) 


Die gesamte Entropieströmung, einschließlich der durch Kon- 
vektion, ist also, wie zu erwarten, gleich 


(26) 10,5, + &,D,. 4 


Auch hier ist zu fordern, daß das Integral der lokalen Entropie- 
erzeugung ı%, über das abgeschlossene System integriert, nicht negativ 
ist. Darüber hinaus verlangen wir jedoch, daß + selbst nicht negativ 
ist. Bei den linearen phänomenologischen Ansätzen des nächsteı 
Abschnitts ist dies schon eine Folge der ersten Forderung, also 
des 2. Hauptsatzes der Thermodynamik. 

Der Ausdruck für die lokale Entropieerzeugung hat eine ein- 
fache Form. Er besteht aus einer Summe von Produkten. Der 
erste Faktor gibt jeweils die strömende Größe an, für die ein Er- 
haltungssatz gilt (Massen-, Energie-, Impulsstrom), der zweite Faktor 
gibt die Ursache für die entsprechende Strömung (äußere K Krak baw. 
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Gradient des Planckschen charakteristischen Potentials, Temperatur- 

gefälle, Geschwindigkeitsgefälle. Mit Onsager werde 


(27) 


> > 
(28) Lp, = K, + T grad ®, 
als „Kraft“, die den Energiestrom und 


29) 1 (5 


ik 2 Ox, Ox; 


als „Kraft“, die den Impulsstrom (innere Reibung) 
Dann ist 


3 3 
(30) T4 = + (D,,L,,) + (Ww, Ly) + a > 
i=1 k=1 
Diesen Ausdruck kann man als Energiedissipation in der 
Volumeneinheit bezeichnen. Das Reibungsglied in (30) ist ja nichts 
anderes als die Rayleighsche Dissipationsfunktion der Hydro- und _ 
Aerodynamik. 
Bei anderen irreversiblen Prozessen ist derselbe Weg gangbar. 
Die „Kräfte“ sind immer als Koeffizienten der entsprechenden 
Ströme in dem Ausdruck TF zu definieren. Das Produkt aus 
„Kraft“ und zugehörigem Strom hat immer die Dimension einer | 
Energiedichte. 


4. Der phänomenologische Ansatz 


sind nur die Erhaltungssätze und allgemeine Beziehungen der Thermo- 
dynamik benutzt worden. Um Probleme der Thermodiffusion durch- 
zurechnen, müssen Annahmen über den Zusammenhang zwischen 
Strömen und „Kräften“ gemacht werden. Den einfachsten phäno- — 
menologischen Ansatz, der den wirklichen Verhältnissen in einem 
weiten Bereich entspricht, stellt eine lineare Beziehung zwischen 
Strömen und Kräften dar. Von der inneren Reibung soll zunächst 
abgesehen werden. Ihre allgemeinste Form ist bei Berücksichtigung 
der Isotropie des Gases 


> > > 
= @,Lo,+ Lp, + 

> > > 
= Go Lp, + Lp, + G33 Ly, 


= dy, Lp, + @, Lp, + Ly. 
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Daß man in diesem Ansatz z. B. einen Wärmestrom mit einem 
Konzentrationsgradienten bzw. einem Diffusionsstrom koppelt und 
ihn nicht nur proportional zum Temperaturgefälle ansetzt, hat zu- 
nächst denselben Grund, aus dem man in einem anisotropen Kristall 
auch eine Komponente des Wärmestroms in der y-Richtung ansetzt, 
wenn der Temperaturgradient in der z-Richtung liegt. Eine tiefere 
Begründung liefert die kinetische Auffassung dieser Vorgänge, wo- 
nach ein diffundierendes Molekül Energie mit sich führen kann. 

Der Ansatz (31) enthält neun Koeffizienten «,,; diese Zahl 
reduziert sich auf höchstens sechs unabhängige Koeffizienten, da 
nach der Definition der Diffusionsströme bei beliebiger Druck-, 


> > 
Konzentrations- und Temperaturverteilung D, = — D, gilt. Daher ist 
(32) = — i=1,2,3. 
Eine weitere Reduktion ergibt sich daraus, daß die Energie- 
> > 
dissipation T+ bei beliebigen Lp; und Ly nicht negativ sein darf. 
Durch Einsetzen von (31) in (30) ergibt sich daher 
> > => >) 
lau Lp, + 02 Lp,; Lp, L»,) + (&ı3 + @3ı) (Lp,, Ly) 
> >> 
+ (032 — 613) (Lp,, Ly) + Ly = 0. 


> 
Damit dieser Ausdruck z.B. für Ly= 0 nicht negativ wird, 
. > . 
muß = — sein; damit er für Ly,= Lp, nicht negativ 
wird, muß 2 =— @;, sein. Dann lautet die Bedingung 7#=0, 
iy (Li, = L»,) + + 31) (Lp, Lp,, Ly + Liv = 0. 
Hieraus folgt noch die Ungleichung 


die einer von Boltzmann’) fiir die thermoelektrischen Effekte 
ebenfalls auf Grund thermodynamisch-phänomenologischer Betrach- 
tungen aufgestellten Beziehung entspricht. 

Die Nützlichkeit des Onsagerschen Reziprozitätssatzes besteht 
darin, daß er die Zahl der unabhängigen Koeffizienten in einem 
solchen Ansatz noch weiter einschränkt; es muß nämlich gelten 


(34) @,= @,, für alle 1, k. 
Damit ist die Zahl der unabhängigen Koeffizienten auf drei 


reduziert. Sie lassen sich durch die Wärmeleitfähigkeit im stationären 


1) L. Boltzmann, Sitzber. d. Akad. d. Wissensch., Wien, Math.-naturw. 
Klasse 96. S. 1258. 1888. 
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Zustand mit DB « - 0, den gewöhnlichen Diffusionskoeffizienten und 
den Koeffizienten Thermodifiusion ausdrücken. 

Der allgemeinste lineare Ansatz im Falle einer Anisotropie, 
wie er bei der Diffusion mit Wärmeleitung in Kristallen zu machen 


wäre, würde in der Darstellung jeder Komponente von D, 3 D,, W durch 


eine lineare Funktion der Komponenten von E; Bi: | bestehen; 
d.h. jeder Koeffizient @,, in (31) wäre durch einen Tensor mit 
9 Komponenten zu ersetzen. Die Beziehung (31) würde dann ins- 
gesamt 81 Koeffizienten enthalten, von denen aber bei Berücksichtigung 
der Symmetrieverhältnisse im Kristall und des Onsagerschen Rezi- 
prozitätssatzes nur ein Teil unabhängig ist. 

Für die Thermodiffusion von paramagnetischen Gasen im Magnet- 
feld ist der Onsagersche Reziprozitätssatz in einer etwas ver- 
allgemeinerten Form anwendbar. Dann ist die Isotropie zerstört 
und es ist der allgemeine Ansatz zu machen, der die neun Kom- 


> > 
ponenten von D,, Dy, W als lineare Funktionen der neun Kom- 
ponenten von darstellt. Das Koeffizientenschema sei 
> 
A: wv =1,2...,9. Dann gilt nach Onsager im Magnetfeld-H 
(H) = 4, (— 
u» \ 


Es läßt sich zeigen, daß es in diesem Fall neun unabhängige 
Koeffizienten gibt, die nur Funktionen des Betrages aber nicht der 


Richtung von H sind. — 

Um den Ansatz (31) mit den Ergebnissen von Enskog zu ver- 
gleichen und um den Zusammenhang der Koeffizienten «,, mit den 
Koeffizienten der Diffusion, Thermodiffusion und Wärmeleitung zu 


erkennen, ist es notwendig, die Größen Fu, =: und i aus (27), 
28) und (29) in (31) einzusetzen. Der Ansatz (31) stimmt dann mit 
den Ergebnissen von Enskog in allen Einzelheiten überein; ins- 
besondere ergibt sich für die Koeffizienten «,,, wenn man sie durch 
die Größen x, A,» der Enskogschen Arbeit ausdrückt: 


7. ae)". “+ RT | 
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Der Onsagersche Reziprozitätssatz (34) wird also durch die 
Enskogsche Rechnung für das Problem der Thermodiffusion he. 
stätigt; insbesondere ist auch &,,= «@,,, da ja bei Enskog, der mit 


einatomigen Molekülen rechnet, u, = 1,2. 


Die Gl. (31) lassen sich unter Berücksichtigung von (32), (34) 
(27), (28), (29) und (18) in folgende Form bringen 


a: [x +Tgrad® —K,—Tg — grad 7 


1 
x” * 


a(u,—U,)—b 
Tz, 


Abkürzungen eingeführt _ 


b = = = — = — 


Der Diffusionsstrom erscheint in fiinf Teile aufgespalten, von 
denen die ersten vier (mit dem Koeffizienten a) von der Verschieden- 


heit der GréBen K, — T grad @,(i = 1,2), also neben den äußeren 
Kräften vom Gradienten der Differenz der Planckschen charak- 
teristischen Potentiale ®, herrühren. Der erste Anteil kommt nur 
bei verschiedenen äußeren Kräften pro Gramm für die beiden 
Molekülarten in Betracht, er fällt also im Schwerefeld und bei der 
Ultrazentrifuge weg. Daß trotzdem eine Trennung von Molekülen 
mit verschiedenem Molekulargewicht in einem solchen Kraftfeld zu- 
standekommen kann, beruht darauf, daß dieses Kraftfeld einen 
Gradienten der Gesamtdichte 9 erzeugt. Dann wird der zweite Anteil 
des Diffusionsstroms wirksam; die schwerere Komponente diffundiert 
nach Stellen höherer Gesamtdichte 9 wegen a > 0 [vgl. (41). Der 
dritte Anteil ist die gewöhnliche Diffusion infolge örtlich ver- 
schiedener Konzentration o,/o. Er wirkt im Sinn eines Konzen- 
trationsausgleichs wegen a>0. Der vierte Anteil liefert einen 
Beitrag zur Thermodiffusion, d. h. zur Diffusion durch das Tem- 
peraturgefälle, falls u, + u, ist; die Moleküle mit der größeren 
inneren Energie pro Gramm wandern nach den Orten höherer 
Temperatur. Einen weiteren Beitrag zur Thermodiffusion liefert der 
fünfte Anteil des Diffusionsstroms in (35), der die weitere phäno- 
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4 a Konstante b enthält. Es sei darauf hingewiesen, daß 
auch der zweite Anteil des Diffusionsstroms einen indirekten Beitrag 
zur Thermodiffusion liefern kann, nämlich z.B. dann, wenn im ganzen 
Gas der Druck konstant ist, aber die Temperatur und damit die Dichte 
örtlich verschieden sind. Hier erzeugt das Temperaturgefille erst 
einen Dichteanstieg; in diesem diffundieren die schwereren Moleküle 
nach Stellen höherer Gesamtdichte, also niedrigerer Temperatur. 

Als Beispiel sei kurz das isotherme atmosphirische Gleichgewicht 


einer Mischung von zwei Gasen besprochen. Hier ist X, =K, = Erd- 
beschleunigung mit dem Betrag g. Der zweite Anteil des Diffusions- 
stroms hebt sich, wie man leicht nachrechnet, gegen den dritten weg, 


da bei diesem Gleichgewicht In —. RT 


Erdboden, o,, = Dichte der Komponente 7 an der Erdoberfläche). Da 


(2 = Hohe über dem 


überdies grad T = 0 angenommen wurde, so folgt D, = 0, wie es im 
Gleichgewicht sein muB. 

Der Diffusionsstrom (35) läßt sich noch in anderer Form dar- 
stellen, indem die Dichte o mittels (16) durch den Gesamtdruck und 
und die Konzentration der Komponente 1 ausgedrückt wird. Es 
ergibt sich 


Hier erscheint der Diffusionsstrom in anderer Weise zerlegt als 
in (35). Diese Darstellung ist dann zweckmäßig, wenn der Gesamt- 
druck p in allen Teilen der Gasmischung derselbe ist, wie es bei 
Diffusionsversuchen gewöhnlich der Fall ist. Den Anteil des 
Diffusionsstroms, der proportional zum Druckgradienten ist, be- 
zeichnet man als Druckdiffusion’). Der Koeffizient des zum Tem- 
peraturgradienten proportionalen Gliedes ist der Thermodiffusions- 
koeffizient bei konstantem Druck. Sein Verhältnis zum Koeffizienten 
der gewöhnlichen Diffusion infolge Konzentrationsgefälles bei kon- 
stantem Druck ist im wesentlichen der Trennfaktor des Thermo- 
diffusionsprozesses. 

Ks ist hier noch auf folgendes hinzuweisen. Die Diffusions-- 
ströme D, (i = 1,2) beziehen sich, wie schon erwihnt, auf eine mit 


1) S. — FERNEN Magazine [6] 34. S. 146. 1917. 
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der Geschwindigkeit v, also mit der Schwerpunktsgeschwindigkeit 
des Volumenelements, mitbewegte Fläche. Wird in einem ab- 
geschlossenen Gefäß durch Herstellung eines Temperaturgradienten 
eine Thermodiffusion bei konstantem Druck eingeleitet, so difiun- 
dieren durch eine im Raum feste Fläche gleichviel Moleküle der 
einen und der anderen Art in entgegengesetzten Richtungen. 
D.h. die mittlere Geschwindigkeit der Moleküle ist hier gleich Null, 
nicht aber die Schwerpunktsgeschwindigkeit, weil sich in einem Teil 
des Gefäßes die leichteren, in einem anderen Teil die schwereren 
Moleküle ansammeln. Es ist daher eine Umrechnung des Diffusions- 
stromes bezüglich der mit der Geschwindigkeit v bewegten Fläche auf 
eine im Raum feste Fläche durchzuführen. Sie liefert das Ergebnis, 
daß bei annähernd gleichem Molekulargewicht der beiden Komponenten 


der relative Unterschied der Diffusionsströme bei mitbewegter bzw. 
(M, — M,) 
(M,+M,) 
Isotopentrennung mit M, = M, in 2 Regel vernachlässigt werden. 


Der Wärmestrom (36) ist ebenfalls zusammengesetzt. Der erste 
Teil ist der Energietransport infolge Diffusion, der zweite Teil ist 
der Wärmestrom infolge des Temperaturgefälles. Bei einer Messung 
der Wärmeleitfähigkeit von Gasmischungen ist dies zu beachten. 


ruhender Fläche von der Ordnung ist. Er kann also bei der 


> > 
Der Einfachheit halber sei vorübergehend K, = K, vorausgesetzt, 
Nimmt man die Messung des Wärmestroms am gleichmäßig durch- 
mischten Gas bei konstantem Druck vor, so ist zu Beginn eine 
Thermodiffusion vorhanden, aber grad (*) =0. Der Wärmestrom 


berechnet sich durch Elimination von D, aus (38) und (36). Er ent- 


man dagegen, bis sich ein stationärer Zustand eingestellt hat, sc 
/ bh? 
T 

und dies ist nach (41) stets positiv. Durch Vergleich dieser beiden 
Wärmeleitfähigkeiten ist es im Prinzip möglich, den Koeffizienten 
der Thermodiffusion zu bestimmen. ‘Für Gasmischungen mit wenig 
verschiedenen Molekulargewichten und gleichem StoBquerschnitt 
(d.h. Kraftgesetz) der beiden Komponenten ist der relative Unter. 
schied der beiden Wärmeleitfähigkeiten von der Größenordnung 
(= - M, , 

‚M,+M, 
Unterschied wesentlich größer werden. 


spricht einer Wärmeleitfähigkeit 


ergibt sich wegen D, = 0 die Wärmeleitfähigkeit zu 


Bei verschiedenen Kraftgesetzen kann dieset 
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Nun noch einige Bemerkungen zur inneren Reibung. Der 
übliche und den Erfahrungen weitgehend entsprechende phänomeno- 
logische Ansatz ist 


8 8 3 
öv öv, 


i=1 m=1 fet 

Die Zahl 81 der Konstanten reduziert sich wesentlich wegen der 
Isotropie; ferner ist wegen des Schwerpunkts- und des Drehimpuls- 
satzes p,, = P,;; der Reziprozititssatz von Onsager fordert 
Gym = Lumix) das gibt aber hier keine weitere Beschränkung. Es 
bleiben zwei unabhängige Koeffizienten, die gewöhnliche Reibungs- 
konstante 7 und die Volumenviskosität; die letztere wird gewöhnlich 
(auch bei Enskog) dadurch eliminiert, daß man den Druck in der 
Gleichung des idealen Gases dem Mittelwert der drei Normal- 
drucke p,,; Pass Pag gleichgesetzt. 

Für die Energiedissipation (30) ergibt sich mit den Ansätzen (31) 
und (39), wenn 4 die Reibungskonstante ist 


> > \2 > > > >? 
Toe =@4 (Lp, + 2b (Lp, Lp,; Ly) _- 
EM 2 + Ov, — $2)" 
dv, Ov,\? , (Or, \*? Ov, 
| + (50, + aa) + (ae, + on) + + |: 
Der Koeffizient von 7, die Rayleighsche Dissipationsfunktion, 
sei im folgenden mit % abgekürzt. Damit das Integral über # 


für ein abgeschlossenes System bei beliebigen Bedingungen positiv 
ist, ist notwendig und hinreichend, daß 


(441) 7>0, a>0, ac—b?>0 und damit auch c>0. 


Also ist «# selbst immer positiv. 
Der Ausdruck (40) läßt eine bemerkenswerte Umformung zu. 


(40) 


> 
Elimination von Bee Lp, mit Hilfe der ersten Gl. (31) gibt 
> 2 

Damit scheint die Energiedissipation additiv aus drei Gliedern 
zusammengesetzt, die auf je einen der irreversiblen Prozesse, nämlich 
Diffusion, Wärmeleitung und innere Reibung, Bezug haben. Es ist 
so, als ob jeder dieser Prozesse seinen spezifischen Beitrag zur 
Energiedissipation gibt und als ob ihre gegenseitige Kopplung, die 
sich ja in der Thermodiffusion äußert, keine gemischten Glieder 
in T # liefert. Eine Verallgemeinerung dieses Ergebnisses ist jedoch 
nicht erlaubt und führt in anderen Fällen sogar zu falschen 
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Folgerungen‘). Das ist auch daran zu erkennen, daß der Aus- 
druck (42) in seiner Zusammensetzung heterogen ist; er enthält teils 


Ströme (D,) , teils „Kräfte“ (grad T) und es besteht natürlich eine 
Willkür, für welche irreversible Erscheinung man den Strom, bzw. 
die „Kraft“ in die Energiedissipation einsetzt. Sind nur zwei irre. 
versible Prozesse in isotropem Material einander überlagert, so ist 
eine Darstellung von «#* ohne Kopplungsglied immer möglich, indem 
für den einen der beiden Prozesse der Strom, für den anderen die 
„Kraft“ eingeführt wird; bei mehr als zwei überlagerten irreversiblen 
Prozessen geht das im allgemeinen nicht mehr; aber auch schon 
bei zwei irreversiblen Prozessen ist eine solche Zerlegung nicht mehr 
möglich, wenn sie unter dem Einfluß eines Magnetfeldes, d. h. unter 
einer Einwirkung rotatorischen Charakters ablaufen). 

Die obigen phänomenologischen Ansätze scheinen zunächst nicht 
die allgemeinsten möglichen zu sein, selbst bei Forderung einer 
linearen Beziehung zwischen Strömen und „Kräften“. Der all- 
gemeinste Ansatz müßte eine lineare Beziehung zwischen Diffusions-, 


Energie- und Impulsstrom einerseits und den „Kräften“ | i 
und L;, andererseits unter Berücksichtigung des Onsagerschen 
Reziprozitätssatzes herstellen. Danach würden die Diffusionsströme 
und der Wärmestrom auch von den Geschwindigkeitsgradienten und 
die Reibungsdrucke auch von den Gradienten der Konzentrationen 


und der Temperatur und von Bi abhängen. Nach Enskog 
besteht eine solche Abhängigkeit bei einatomigen Gasen nicht. Es 
läßt sich aber auch allgemein auf Grund der Isotropie der Gas- 
mischung und wegen p,, = p,, einsehen, daß ein solcher Zusammen- 
hang bei Mischungen beliebiger idealer Gase nicht möglich ist. 


5. Thermodiffusion in realen Gasen und Flüssigkeiten 


Die Überlegungen zur Thermodiffusion mit Konvektion in 
Mischungen idealer Gase lassen sich ohne weiteres auf Mischungen 
realer Gase und auf flüssige Lösungen "mit zwei oder mehr un- 
abhängigen Komponenten übertragen. Die Kontinuitätsgleichungen 
für die Mengen der Komponenten, für den Impuls und für die Energie 
bleiben in der in (7), (8), (9) bzw. (12) angegebenen Form bestehen. 
Die Energie U und die Entropie S lassen sich aber, außer bei ver- 
dünnten Lösungen und realen Gasen mit geringer Konzentration 
einer Komponente, nicht mehr in der Form (13) und (14) darstellen. 
Sie sind in allgemeiner Form durch die Zustandsgrößen auszudrücken. 
Als solche seien neben der Temperatur T die Dichten o,@ = 1,2,...," 


1) J. ere 2.2.0. 
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der einzelnen Komponenten DENN, n sei die Zahl der Komponenten. 
Für jede ist ein Diffusionsstrom D. einzuführen und es gilt 


(43) D+ Dy 


n 


(44) TdS=dU — 


i=l 


Sie sind auf 1 g bezogen, ebenso wie die thermodynamischen 
ale u, der i. die 0; 


gegen Druck p Also ist 


von U und o, nach der Zeit aus den Gl.(7) und (12) eingesetzt 


werden (in (7) ist t= 1,2, ..., 
= DE, B) + Sp, aivd 
5) 
(45) 


Diese Gleichung hat die Form einer Kontinuitätsgleichung vom 
Typ (22), wenn 


(46) 


als Entropieströmung (ohne RR und % 


| 
| Th = > (D, T grad | grad T 


(47) 
Ov; 
| > (P9;,— Pit) 


kal 


als Energiedissipation definiert wird. 

Damit ist gezeigt, daß die Ergebnisse (23) und (24) in diesem 
allgemeinen Fall bestehen bleiben. Auch hier ist der phänomeno- 
logische Ansatz (31) und (39) unter Berücksichtigung des Onsager- 
schen Reziprozitätssatzes zu machen und alle Folgerungen der letzten 
Abschnitte, die nicht auf der Anwendung der Gesetze der idealen _ 
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(Gase beruhen, behalten ihre Gültigkeit. Bei zwei Komponenten ist 
die Zahl der unabhängigen Koeffizienten dieselbe wie im Fall der 
Mischung zweier idealer Gase. 

In festen Lösungen und Mischungen (Legierungen) ist eine 
eigentliche Konvektion auszuschließen. Man wird hier die Diffusions- 
ströme nicht auf eine mit der Geschwindigkeit © bewegte Fläche, 
sondern auf eine im Gitter feste Fläche beziehen. Dann ist aber 


n 
nicht mehr I" D,= 0. Statt der Konvektion sind die Verschiebunger 
i=1 
der Gitterpunkte, statt der Reibungsdrucke die elastischen Spannungen 
einzuführen. Die phänomenologischen Ansätze für die Beziehungen 
zwischen den Strömen und den „Kräften“ (einschließlich der inneren 
Spannungen) sind, solange man mit linearen Ansätzen durchkommt, 
im wesentlichen dieselben wie oben; auch der Onsagersche Rezi- 
prozitätssatz ist wieder anwendbar. Die Zahl der unabhängigen 
Koeffizienten ist jedoch wegen der Anisotropie (außer bei kubischen 
Kristallstrukturen) größer als bei Flüssigkeiten und Gasen. 
. Zusammenhang mit den Überführungswärmen 
System sei stationär genannt, wenn es 
mit seiner Umgebung nur Wärme, aber keine Stoffmengen aus- 
tauscht, und wenn sich in seinem Inneren keine zeitlichen Ände- 
rungen der Zusammensetzung und der Temperatur abspielen. Di 
Temperatur braucht jedoch nicht überall dieselbe zu sein. Zur Be- 
handlung solcher stationärer Systeme hat Eastman’) den Begriff 
der Überführungswärmen eingeführt. Wagner?) hat in sehr all- 
gemeinen Untersuchungen diesen Begriff insbesondere auch auf 
stationäre Systeme mit elektrischen Erscheinungen angewandt. Beim 
Problem der Thermodiffusion ist die -Definition der Überführungs- 
wärme so auszusprechen. Es seien zwei gleich zusammengesetzte 
Mischungen von zwei Gasen mit gleicher Temperatur gegeben. Nun 
soll 1g des Gases 1 vom ersten in das zweite System, 1 g des 
Gases 2 vom zweiten in «as erste System quasistatisch übergeführt 
werden (die beiden Mischungen sollen so große Massen haben, dab 
sich bei der Überführung von je 1 g die Konzentrationen praktisch 
nicht ändern. Um den Vorgang isotherm durchzuführen, ist eine 
Überführungswärme Q dem System 1 zu entziehen und dem System 2 
zuzuführen. In unserer Darstellung des Diffusionsvorganges ist also 


1) E.D. Eastman, Journ. Amer. Chem. Soc. 48. 8. 1482. 1926; 50. S. 28 
292. 1928. 
2) C. Wagner, Ann. d. Phys. [5] 3. S. 629. 1929; 6. S. 370. 1930. 
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die Überführungswärme nichts anderes als die Wärmemenge, die 


-> > ‘ 
mit den Diffusionsstrémen D, und D, mitstrémt, wenn insgesamt 
je 1g der beiden Komponenten in entgegengesetzten Richtungen 


> 
durch eine mit der Geschwindigkeit » bewegte und zu ihr senk- 
rechte Fläche diffundiert ist. Nach (36) ist für isotherme Diffusion 


> 
a 


W = 
Die Überführungswärme Q ist daher durch den Koeffizienten b/a 
gegeben. Die Bedeutung der Überführungswärme liegt darin, daß 
das stationäre Gleichgewicht (zu unterscheiden vom thermodyna- 
mischen Gleichgewicht mit einheitlicher Temperatur des Systems) 
durch sie ausgedrückt werden kann. Im stationären System ist 
D,=D,=0. 
Aus (35) folgt daher 
(48) K, — K, + Tgrad(®, — 7 grad | 
Für zwei Gasmischungen, die durch eine Röhre verbunden sind ce = 
und die durch Wärmebäder auf den Temperaturen T und T + dT ge- Sr é 
> 
halten werden, stellt sich im Falle konservativer Kräfte (A, = — grad V ) en ur 
ein Konzentrations- und Druckunterschied ein derart, dab E 


(49) —V,)+ Td@, — ®, = “aT. 
Diese Gleichung reicht gerade aus, um den Konzentrations- 
unterschied bei gegebenen dT zu bestimmen, denn der Druckunter- 


= 
u 


schied ergibt sich aus (8). Für das stationäre System ist ja v=0 
ud -dp=o,dV, +0,dV,. 

Daß die Überführungswärme das stationäre Gleichgewicht be- 
stimmt, beruht in unserer Behandlung des irreversiblen Vorgangs 
der Diffusion und Wärmeleitung darauf, daB «,=«,, [vgl. (31) _ 
und (34)] und damit auf dem Onsagerschen Reziprozitätssatz. 

Eastman und Wagner leiten Beziehungen vom Typ (49) mit 
Hilfe von quasithermodynamischen Betrachtungen ab. Bei Eastman 
ist speziell die Annahme gemacht, daß der Entropieaustauch 
zwischen nichtisothermen Teilsystemen nur durch Stoffübergang, 
aber nicht durch Wärmeleitung erfolgt. Bei Wagner ist dieselbe — 
Annahme in einer anderen Form gemacht: er nimmt an, daß der 
Austausch der Komponenten zwischen zwei Teilsystemen bei einem 
gedachten Kreisprozeß in so kurzen Zeiten erfolgt, daß die durch 
Diffusion übertragene Wärmemenge groß gegen die durch Wärme- 
leitung im stationären Zustand übertragene Wärmemenge ist. 

Annalen der Physik. 5. Folge. 39. 24 
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Daß diese Annahme im vorliegenden Fall auf die Beziehung 
. . . . . 
&;=@;,, führt, ist leicht einzusehen. FürD, = D, = 0 ist ja nach(31 
und (32) und wegen der mit rein thermodynamischen Betrachtungen 
abgeleiteten Beziehungen @,, =— a@,, >0, a4, =— @,, (vgl. S. 344) 
Oso & =o 
W = (a, “2 ) Ly. 


Vernachlässigt man die stationäre Wärmeleitung, so setzt man damit 
233%, = %, %, und das ist mit(33) nur verträglich, wenn «,, = «,, ist, 

Zur Ableitung dieses Ergebnisses ist jedoch keine quasithermo- 
dynamische Überlegung erforderlich, es folgt ja unmittelbar aus dem 
Onsagerschen Reziprozitätssatz. Überdies ist diese Überlegung 
mit der Annahme vernachlässigbarer stationärer Warmeleitfihigkeit 
nicht ohne weiteres gerechtfertigt, weil Wärmeleitung und Diffusion 
nicht unabhängige irreversible Prozesse sind; im Gegenteil sind ihre 
Mechanismen eng miteinander verknüpft. Wagner hat die Schwierig- 
keiten einer solchen Annahme ausführlich erörtert (a. a. O. S. 643), 
Daß solche Überlegungen trotzdem häufig zu richtigen Ergebnissen 
führen, bedeutet natürlich für sie eine gewisse Rechtfertigung. Da 
sie in manchen Fällen aber auch versagen können [z. B. bei der 
analogen Behandlung der thermoelektrischen Effekte im Magnet- 
feld})], ist es zweckmäßiger, sie durch Anwendung des weiter reichen- 
den Reziprozitätssatzes von Onsager zu ersetzen. 

Auch bei einer Mischung von beliebig vielen idealen oder real 
Gasen oder flüssigen Lösungen (die Zahl der Komponenten sei 
reichen die Überführungswärmen aus, um das stationäre Gleichgewicht 
festzulegen. Für die Diffusionsströme gilt (43). Der phänomenologische 
Ansatz (31) ist zu verallgemeinern in 


n 


5 
(50) = > lox P;Lw; 
k=1 


n 

(51) > bilo: +7Ly. 

k=1 

Die 5 sind wie in (28) bzw. (47) [vgl. definiert. 
Onsager gilt 


(52) = Bi; 
ferner gilt wegen (43) 


(53) > = 0 für k= 1,2,..., 


i=1 
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Q, sei nun die überführte Wärme, wenn 1 g der Komponente k 
(k=1,2,...n) in der einen Richtung, 1 g der Komponente n in 
der entgegengesetzten Richtung überführt wird; per def. ist also 


Q,=0. Wir behaupten, dab 3; = ze @,,Q,. Wird dies nämlich 
%=1 
in (51) eingesetzt und werden dann die Lp; mittels (50) ausgedrückt, 
so entsteht wegen (52) 


n n n 
> —> y | 
(4) > Q,D, + (7 > > 0,0, Q, Lw. 
k=1 i=l kel 
Andererseits gilt für den stationären Zustand (D, = 0), wenn 
man die 9, in (50) durch die Überführungswärmen ersetzt a 
é em 2 IHR, 
Y 
k=1 
Diese n homogenen Gleichungen sind nicht voneinander un- 


abhängig. Man kann also nicht auf om + QO: Re = 0 schließen. 
Wegen (53) kann man sie jedoch in der Form schreiben 


n—1 
Q, Lw) 0; i= 1, 2,. — 1. 


k= 1 


Diese n —1 homogenen Gleichungen für die n— 1 Funktionen 


Lp, — Ly; + Q, Lw sind voneinander unabhängig; es folgt also 
re > > 
— Loa + Q, Lw = 9; km 1,2,... n—1. 


Diese Bedingungen reichen zur Festlegung des stationären 
Gleichgewichts aus. Sie sind in den von Wagner angegebenen 
allgemeineren Gleichungen enthalten. 

Die Überlegungen dieses Abschnitts lassen sich mit geringen 
Änderungen auch auf feste Lösungen übertragen (vgl. den letzten 

Absatz des Abschnitts 5). 
Zusammenfassung 

Zur Behandlung von irreversiblen Vorgängen ist i in vielen Fällen 
die sogenannte thermodynamisch-phänomenologische Methode geeignet. 
Ihre thermodynamische Seite besteht in der Aufstellung der Er- 
haltungssätze für die Stoffmengen, für Energie und Impuls, bei 
Systemen mit elektrischen Erscheinungen auch für die Ladung und 
in der Aufstellung der Entropiebilanz für jedes Volumenelement. 
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Der zweite Hauptsatz der Thermodynamik ist bei irreversiblen Vor- 
gängen so zu formulieren, daß eine Entropieströmung existiert, die 
räumlich verteilte Quellen (lokale Entropieerzeugung), aber keine 
Senken hat. Die phänomenologische Seite der Methode besteht in 
sogenannten phänomenologischen Ansätzen, die lineare Beziehungen 
zwischen den Strömen (Energie-, Impuls-, Diffusionsstrom usw.) und 
den sie hervorrufenden Ursachen (Temperatur-, Geschwindigkeits-, 
Konzentrationsgradient usw.) vermitteln. Diese Ansätze geben im 
allgemeinen die wirklichen Verhältnisse in einem sehr weiten Bereich 
wieder. Besonders zweckmäßig ist es, für die Ursachen der ver- 
schiedenen Ströme geeignete verallgemeinerte Kräfte einzuführen: 
so ist der Diffusion als verallgemeinerte Kraft nicht der Konzen- 
trationsgradient, sondern im wesentlichen der Gradient der Planck- 
schen charakteristischen Potentiale zuzuordnen. Eine sehr wichtige 
Rolle für diese phänomenologischen Ansätze spielt der Onsagersche 
Reziprozitätssatz, der auf allgemeinen schwankungstheoretischen Über- 
legungen beruht und der für irreversible Vorgänge eine wichtige 
Ergänzung der beiden Hauptsätze darstellt. Er gibt Symmetrie- 
beziehungen für die Kopplung, die bei mehreren einander über 
lagerten irreversiblen Prozessen eintritt. Mit dieser Methode lasseı 
sich nicht nur stationäre, sondern zeitlich beliebig ablaufende Vor- 
gänge behandeln. Zunächst wird sie am Beispiel der Wärmeleitung 
in anisotropen Kristallen auseinandergesetzt. Als Anwendung wird 
ausführlich die nichtstationäre Thermodiffusion mit Konvektion in 
einer Mischung von zwei idealen Gasen behandelt. Auf die Ver- 
hältnisse bei der Thermodiffusion von paramagnetischen Gasen im 
Magnetfeld wird kurz hingewiesen. Die Überlegungen lassen sich 
leicht auf die Thermodiffusion in realen Gasen, Flüssigkeiten und 
festen Körpern ausdehnen. In stationären Systemen (ohne Diffusions- 
ströme, mit zeitlich konstanter Temperaturverteilung) stimmen die 
Ergebnisse mit den von Eastman und Wagner mit Hilfe des 
Begriffs der Überführungswärme auf quasithermodynamischer Grund- 
lage gewonnenen Ergebnissen überein; diese Überführungswärmen 
haben im phänomenologischen Ansatz eine einfache und anschauliche 
Bedeutung. Gewissen Einwänden, denen die quasithermodynamische 
Methode ausgesetzt ist, unterliegt die thermodynamisch-phänomeno- 
logische Methode nicht. 


Berlin, I. Institut für theoretische Physik und I. Physikal ikalisches 
Institut der Universität Berlin. 


(Eingegangen 7. April 1941) 
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in 
Von Karl Bechert 


Eine Zusammenstellung von Bemerkungen, die fiir die us: von 
Problemen der Wellenausbreitung in Gasen nützlich sein können. — § 1. Ein- 
leitung, Ergebnisse. $ 2. Über die Ausnützung von Ähnlichkeitstransformationen. 
$3. Die Veränderlichen m, ¢. § 4. Eine nicht-lineare partielle Differential- 
gleichung zweiter Ordnung und ihre Integration durch „Aufspaltung“. § 5. Über 
das Problem der Ausbreitung ebener Wellen in idealen Gasen, wenn der 
Zusammenhang polytroper Zustandsänderung nicht vorausgesetzt wird. 


$1. Einleitung, Ergebnisse 


polytroper Zustandsänderungen vollständig erledigen, 

Riemann gezeigt hat; für die physikalisch meist vorkommenden 
Fälle kann die Lösung sogar in handliche Form gebracht werden’). — 
Für die Gleichungen der Ausbreitung von Zylinder- und Kugel- 


Die Behandlung solcher Aufgaben führt auf Fragen und Ergebnisse, a 
die gegenüber den sonst in der theoretischen Physik vorkommenden = 
Rand- und Anfangswertaufgaben partieller Differentialgleichungen = 
ungewohnt sind. Zum Beispiel ist bei ebenen Wellen eine a = 
kehrungsrechnung“ zur Lösung der Anfangswertaufgabe nötig (I, § 6): _ 

für ein Verhältnis x der spezifischen Wärmen: = 3 hat man 

zwei Funktionen F und G zweier Veränderlichen £, 7 aus —. 


Gleichungen zu bestimmen [I, (51)]: 


v 

/ 3 0 
(la) Vz (= | ( 2 ( 
(1b) 0= (5; F(E ++ G(E— 
& 


1) K. Bechert, Ann. d. Phys. [5] 37. 5. 89. 1940; Erweiterung in Ann. eae: 
d. Phys. [5] 38. $. 1. 1940. Im folgenden als Arbeit I und i zitiert. as ce 

IK. Bechert, Ann. d. Phys. [5] 39. S. 169. 1941; h 
Arbeit III zitiert. 
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wobei die Größen E, i physikalisch die „Anfangswerte“ von 4 
darstellen und vorgegebene Funktionen von x sind. Die Lösung dieger 


Aufgabe fiir <=. und die des entsprechenden kompliziertere 


‘ 


Problems fir x = — — dort treten zweite Differentialquotiente 


auf und in verwickelterer Form als die ersten in (1) — ist in d 
Arbeit I gegeben; es dürfte ganz interessant sein, ist aber nicht die 
Absicht der vorliegenden Arbeit, den allgemeinen Algorithmus auf- 
zufinden, der solche Aufgaben in einfacher Weise löst. 

Bemerkenswert scheint mir auch, daß stehende ebene Wellen 
endlicher Amplitude zwischen festen Wänden nicht durch den nahe- 
liegenden Ansatz f(z)g(t) für die Dichteschwankung 9 — o, gegen- 
über dem ungestörten Wert 9, oder die Strömungsgeschwindigkeit u ge- 
liefert werden. Die einzige physikalisch zu stellende Bedingung ist u= 
am Rand, und diese ist mit dem genannten Ansatz hier nicht veı 
träglich. Trotzdem hat u Knoten, also Stellen, an denen es fi 
alle ¢ Null ist; diese Stellen liegen übrigens genau da, wo sie nac 
der üblichen Rechnung, die beliebig kleine Amplituden voraussetz 
liegen sollen. Auch daß die Extrema von o an den Knotenstelle 
von % liegen, ist für endliche Amplituden noch richtig (I, § 7, 2 
dagegen hat o—o, nach der Rechnung für endliche Amplitude 
keine Knoten, im Gegensatz zur üblichen Theorie beliebig kleine 
Amplituden. 

Ergebnisse. Ich will hier von Methoden sprechen, die zur Auf- 
findung spezieller Lösungen solcher und verwandter hydrodynamischer 
Aufgaben dienen können. In $ 2 ist von den Differentialgleichungen 
die Rede, die in den Arbeiten I—III (über ebene und über Zylinder- 
und Kugelwellen) behandelt wurden. Diese Gleichungen erlauben eine 
einfache Ähnlichkeitstransformation. Man kann sie dazu ausnützen, 
spezielle Integrale zu finden, und z.B. die in der Arbeit III ge- 
gebene Lösung mit ihrer Hilfe ableiten ($ 2). Die Methode ist ver- 
allgemeinerungsfähig. 

§ 3 will darauf aufmerksam machen, daß es zweckmäßig sei 
kann, z.B. bei ebenen Wellen statt der Veränderlichen x,t d 


Größen m, t einzuführen; dabei bedeutet m = [ode die in ein 


Zo 
bestimmten Gasschicht enthaltene Masse. Diese Substitution verdanke 
ich Herrn Dr. H. Marx. 
Die Differentialgleichungen fiir die ebene Bewegung lassen sic 
auf eine nicht-lineare Differentialgleichung zweiter Ordnung fiir ew 
Unbekannte zurückführen. Da die Lösung für ebene Wellen bi 
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angenehm aussehende Differentialgleichung streng integriert werden 
kann. Die Art des Rechenverfahrens, das dabei angewandt werden 
muß, ist recht merkwürdig: es kommt darauf hinaus, die Gleichung 
in zwei niedrigerer Ordnung aufzuspalten und diese zu integrieren. 
Das Integral enthält zwei willkürliche Funktionen, obwohl es eine 
Beziehung zwischen ersten Differentialquotienten ist ($ 4). Die Au- — 


gleichung mit zwei unabhängigen Veränderlichen gibt ein erstes 
Integral der Wellengleichung, das zwar auch aus der allgemeinen 
d’Alembertschen Lösung leicht abgeleitet werden kann, das aber 
doch überrascht, wenn man es zum erstenmal sieht. Dabei ergibt 
sich die Frage, wieviel willkürliche Funktionen in einem Integral 
von partiellen Differentialgleichungen vorkommen müssen, damit es 
als allgemeines Integral betrachtet werden kann. Diese Frage ist 
grundsätzlich an Hand der Taylorschen Reihe beantwortbar (§ 4). 
Die Methode der Variablen m, t läßt sich auch dann noch an- _ 
wenden, wenn die Voraussetzung des polytropen Zusammenhangs 
zwischen Druck p und Dichte o nicht von vornherein gemacht wird, 
wenn also nicht vorausgesetzt wird, dab: 
(2) +p, 


n 


sei, wobei a, n und p, Konstante sein sollen. Die allgemein für 
ebene Wellen in idealen Gasen mit Reibung und Wärmeleitung 
geltenden Gleichungen werden in § 5 angeschrieben und auf m,t 
umgerechnet. Für ideale reibungsfreie nicht-wärmeleitende Gase 
ergibt sich als Integral der Gleichungen zwar formal wieder die 
Beziehung (2) mit p,=0, n=x; aber an die Stelle von a tritt 
eine Funktion, die von m in willkürlicher Weise abhängt. Mit 
einem einfachen Eliminationsverfahren gewinnt man eine nicht- 
lineare Differentialgleichung zweiter Ordnung, und mit der in § 2 
beschriebenen Transformationsmethode eine Lösung, welche die 
Zustandsgrößen in Abhängigkeit von einer willkürlich bleibenden 
Funktion von m enthält (85). Eine Diskussion der Lösung zeigt, 
daß die damit beschriebenen Wellen ihre Form beim Weiterwandern 
verzerren — was nach den Ergebnissen der früheren Arbeiten zu er- 
warten war. Theoretisch können also auch hier Stoßwellen auftreten. 


$2. Über die Ausnützung von Ähnlichkeitstransformationen 


Wir gehen von den Differentialgleichungen aus, die für die 
Ausbreitung von ebenen Wellen oder von Zylinder- oder Kugel- 
wellen in reibungsfreien nicht-wärmeleitenden Gasen bei polytroper 
ars, 
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Zustandsiinderung gelten. Es wird also die Gültigkeit von (2) mit 
konstantem a, n, p, vorausgesetzt. Volumkräfte sollen nicht wirken, 
An Stelle von o, wu führen wir wie in der Arbeit III ein: 7 


(3) n=uyn-1; 


k gibt sozusagen die Dimension der Wellenausbreitung an, es ist 
k=1 für ebene Wellen, k = 2 für Zylinderwellen, k = 3 für Kugel- 
wellen. r bedeutet für k = 1 eine kartesische Koordinate (2), fürk = 2 
die Entfernung von der z-Achse eines Zylinderkoordinatensystems, 
für k = 3 die Entfernung vom Koordinatenanfang eines Kugelkoor- 
dinatensystems. Die hydrodynamischen Bewegungsgleichungen werden 
dann durch die eine Gleichung dargestellt: 


die Kontinuitätsgleichung durch: 


é, ‚n-1 . k-1)Vn-1 
&n,+ 


VYn-1 2 


- 


(4b) &+ 


Die Veränderlichen der Gl. (4) kann man einer Ähnlichkeitstransfor- 
mation (Streckung) unterwerfen, ohne daß die Gleichungen sich 
ändern. Setzt man nämlich: 


(da) t=a,t 


mit konstanten Koeffizienten «,, so gelten für die — 
Größen wieder die Gl. (4), wenn nur: 


a, 
as a 


wie das Einsetzen von (5a) in (4) zeigt. ER Gl. (4) sieht man 
außerdem die Transformation t > t-+ const an, die uns aber 
hier nicht weiter beschäftigen soll; physikalisch ist sie übrigens 
selbstverständlich.] Die letzten Bedingungen verlangen nur, dab 


und dab a, = + @, sei. «, und «, bleiben willkürlich; wir 


können also behaupten, daß speziell durch die Transformation: 
(5b) r—Bßr; t > »y& n- 
/ 


die Gl. (4) bei beliebigem konstanten f, y nicht geändert werden. 
Das bedeutet: das allgemeine Integral von (4) muB bei dieser Trans- 
formation formal bleiben. 
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K. Bechert. Ditferentialgleichungen tor 
Wir wollen das an einem einfachen Beispiel überlegen. Die En 

Wellengleichung in zwei unabhängigen Veränderlichen: 

(6a) 

erlaubt die Transformation 


(6b) t=+at, 
mit beliebigem konstanten «, 2. Das allgemeine Integral: 
(6e) 


bleibt in der Tat bei dieser Transformation seiner Form nach er- 
halten; es ist nämlich nach (6b, ce): 


U’ — els(“ * + +N)+g(—t). 


Wir untersuchen jetzt die Folgerungen aus der sehr viel weiter- 
gehenden Forderung, daß die Funktion (6c) bei der ung 
(6b) in sich übergehen soll. Sie bedeutet, daß: 


(Ge) U' t’') = U(2’,t’) = 3 t) 


sein soll, daß nämlich die Funktion U’ aus U einfach dadurch 
hervorgeht, daß man in U die Veränderlichen z,t durch x',t' ersetzt. 
Dann muß zugleich U’ = $ U werden, wie (6e) zeigt. Das erzwingt 
eine Beziehung zwischen @ und /, denn es muß sein (wir begnügen 
uns mit dem oberen Vorzeichen bei f'): 


(7) U’'=U (ea, at) = U(z,1). 


Der Zusammenhang zwischen @ und 3 kann aus (7) durch Differen- 
tiation nach & gefunden werden. Natürlich wäre er auch leicht 
aus (6c) zu finden, aber die Kenntnis des Integrals wollen wir 
jetzt nicht voraussetzen. Es ergibt sich: 

> U, U (a, t); 
und weiter: 


(8) U,-2+U,-t 


Da U’ nur von 2’, t’ abhängt und «, 8 explizit nicht enthält, so 


folgt: — const = !; 8=a'. Es gibt also Lösungen mit der 


Eigenschaft: 
(9) U (az, at) = «!-U(z,t): 
« beliebig. Man kann sie finden, indem man aus (8) auf U = z'- Fie ® 


schließt (die Striche sind jetzt unnötig), und F aus der W ellen- ion 
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u. (6a) bestimmt; dabei ergibt sich eine gewöhnliche Diffe- 
rentialgleichung zweiter Ordnung mit der unabhängigen Veränder- 
lichen x/t. Durch die Forderung (6e) sind eben spezielle Lösungen 
aus der Gesamtheit aller möglichen Lösungen ausgesondert worden, 
In unserem Beispiel ist das allgemeine Integral bekannt, die Methode 
kann daher nichts grundsätzlich Neues über die Integrale geben. 
Wenn es sich aber darum handelt, von einer verzwickt ausselienden 
Differentialgleichung irgendwelche Integrale zu finden, um einen wenn 
auch noch kleinen Einblick in ihre Möglichkeiten zu bekommen, so 

ist die Methode recht brauchbar. 

Zum Beweise leiten wir die in der Arbeit III gefundene Lösung 
aus den Transformationseigenschaften (5b) ab. Wir schreiben: 
(10) n=f(r,t; 

f. 9 sind noch unbekannte Funktionen. Statt der aus (5b) fließenden 
Bedingung, daß die Gl. (10) bei der Transformation (5b) in Glei- 

chungen derselben Art übergehen sollen, stellen wir die schärfere 
Bedingung — sie entspricht (6e) —, daß: 


yf) = Br, t); 79 (rt) =9 (87, 


I 


sei. Damit haben wir auf das allgemeine Integral verzichtet und 
uns auf eine Klasse von Lösungen beschränkt. Wie sie aussehen, 
wollen wir jetzt genauer bestimmen. 

(11) bedeutet einen Zusammenhang zwischen # und y; auf 
dieselbe Weise wie im vorigen Beispiel bekommen wir aus der 


ersten Gl. (11) durch Differentiation nach y: 
(12) fief, (3 1). 
Nochmalige Differentiation nach 7 zeigt, daß notwendigerweise gelter 


muß, entweder (3), +0, also: 


2a) + = 


Gy 7 
4 =const=1; f = -r+(—f, 


= erste Möglichkeit (12a) gibt (die Striche lassen wir jetzt 
der ersten Gleichung: f =F, (4) und aus der zweit 


1 


‚F: (r). Einsetzen der ersten Form in die zweite zeigt, dab 


” .const: das heißt: 


862 ! 
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(13a) n= , const; | 


für die zweite Gl. (11) folgt ebenso: 


(13b) = - const. 


(13) entspricht der singulären Lösung (22a) der Arbeit III, denn t _ 

kann durch ¢+ const ersetzt werden. Zur Bestimmung der Kon- | 

stanten hat man noch (13) in (4) einzusetzen. 
Die zweite Möglichkeit (12 b) gibt: # = const -y! und, wenn wir 


] 
die Striche wieder weglassen und - an Stelle von 1 schreiben: 


(14a) rey =F, (, + const 


Die zweite Gl. (11) hat dieselbe Form wie die erste, daher gilt auch _ 
(hier zeigt sich der Nutzen der Einführung von Veränderlichen, die 
sich in gleicher Weise transformieren: &, nicht 9, transformiert sich 

so wie um n): 


(14b) I 4 (; ; ® 


Daraus folgt erstens, daß r¢& Funktion von r*7 ist — das war der 
allgemeinste Ansatz in III, vgl. dort Gl. (38) — und zweitens, daß 


a+l 


r¢n nur von abhängt — das hatten wir in III durch 


t + const 
Einsetzen in die Differentialgleichungen (4) gefunden, vgl. dort G1. (39), 
(41b), (43a). Auch, daß der Wert a = —1 eine Sonderstellung ein- 
nimmt, ist deutlich; vgl. III, Gl. (43b). 


denn der Rechnungsgang ist derselbe. 
Die Methode ist verallgemeinerungsfähig. 
§ 3. Die Veränderlichen im, ¢ 


Die Differentialgleichungen für ebene Wellen können durch 
Einführung neuer Veränderlicher an Stelle von 2, t a 
gemacht werden. Wir schreiben die Gleichungen in der Form: 


- du 1 d 
(15) pu =0: 
= 9%; dt ou, =0; 


schreiten mit der bewegten Materie. Es sei 


n 
le 
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en 
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nd 
Zur Bestimmung der Funktionen F,, F, hat man (14a,b) in (4 
einzusetzen. Dabei ergibt sich natürlich die in III gefundene Lösung, _ 
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die Masse, die in einer Materieschicht zu bestimmter Zeit vorhanden 
ist (in der y- und z-Richtung sei das betrachtete Volumen je 1 cm 
lang). m hängt von x (und über o auch vont) ab; wir können daher 
an Stelle von z,t auch m,t als unabhängige Veränderliche benützen, 
Auf die Zweckmäßigkeit dieser Substitution hat mich Herr Dr, 
H. Marx aufmerksam gemacht. Wenn wir das Weiterwandern der 
Materieschicht untersuchen wollen, so heißt das, daß wir m dabei 
festhalten, und daß sich x, und x im Laufe der Zeit ändern: 

A dm 

(16 b) | Fr 
andererseits ist bei festem ft: 
(16e) m,=0. 


Die Funktionen, die bei der Variableninderung entstelien, kenn- 
zeichnen wir durch Uberstreichen. Es ist: ; 
F (a,t) = F 
insbesondere: = z(m,t). Aus (16d) folgt: 

(16e) m; 
ferner gilt: 


N IF 
(168) ), = const 


ei denn dF/dt bedeutet seiner Definition nach die Änderung von F 
fiir eine bestimmte Materieschicht. In (16d) kann man links 
_ wieder z = Z(m,t) einsetzen und erhält: 
(16g) F (z(m,t),t) = Pim,t), 
und daraus: 
F, =f. 


ae zweite Gl. (16h) ist die bekannte Beziehung: ae =F u+F, 

0=- = m,u+ m,; damit ist auch m, bekannt: 


denn es ist & = u. Setzt man F = m, so ae wegen (16b): 


m,=— ou. 
- Mit Hilfe von (16c, e,j) kann von den Unabhängigen z,¢ auf die 


Unabhängigen m,t umgerechnet werden. 
Aus den Gl. (15) wird so: 


= 


| 
2 3 
a b 
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| 
a, Aus den ersten Gl. von (16e, h) und aus (16c) folgt: 3 
D 
- 
(f) 
17) = 0; 
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p sei als Funktion von o gegeben, den speziellen Zusammenhang (2) 
brauchen wir nicht vorauszusetzen. Wir nehmen an, die Lösung 
von (17) sei in der Form gefunden: 

(18a) i=g, md: o= g,(m,t). 

Die Umkehrung sei: 

(18b) m=f,(0,@); (t=f,@,@); 

da aber 7, = @ ist, so gilt: 


t 
(18¢) -[yı (m, t)dt = ®(m,t); 


die Integration ist bei konstantem m auszuführen. Durch Einsetzen 
von (18b) in (18c) folgt als Lösung in den physikalisch gegebenen 
Veränderlichen: 
(184) (5,69, 5,6%); t=f,@,a). 
Statt (1Sc) kann man wegen (16i) auch schreiben: 


(18 e) -[ (mn, 8) dm = @(m,t); 


die Integration ist diesmal bei konstantem ¢ auszuführen. 


Zusatz bei der Korrektur am 25. April 1941: 


Liegt die Lösung von (17) unmittelbar in Form der „Umkehrung“ (18b) 
vor, so ist die Auflösung nach x einfacher als in (18c,d). Wenn nämlich m 
und ¢ als Funktionen von u, 9 bekannt sind, so gilt: 


(a) m=m(i,0): t = ¢t(u, 9): 
es ist aber allgemein: ne 

fi m | 


Faßt man andererseits x,t als Funktionen von u,o auf, so folgt: 


t 


A, A, 4,’ 
Durch Vergleich ergibt sich wegen (16e): 
(d) 
und aus (16h): 


(e) 
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wobei das Integral bei festem g auszuführen ist. — Daß sich aus (e) und ( 


das gleiche x,, ergibt, ist eine Folge der Gleichung: sya 
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my 

v @ 


— + 4, = 0, also mit der zweiten Gl. (17). Daher 


die gleichbedeutend ist mit =, 


oO 
ist (g) wirklich das Integral von (e) und (f). — Die rechte Seite von (g) ist nach 
den Voraussetzungen (a) als Funktion von u, 9 bekannt und offenbar ohne 
weitere „Umkehrung“ berechenbar. Man hat damit x und ¢ in Abhängigkeit 
von u, 9. 

Auch beim Vorliegen der Lösung in der Form (18a) kommt man übrigens 
ohne Umkehrung aus, wenn man (18a), (18c) als Darstellung von 7, 7, % durch 
den Parameter m und die Koordinate ¢ auffaßt. 

Zur Erläuterung der Transformation auf m,t: die Gl. (16b) ist gleich- 
bedeutend mit der zweiten Gl. (15), denn beide drücken die Massenerhaltung 
aus. Mathematisch folgt die Gleichwertigkeit aus: 


z 
d 
(i) ie = 0) de J|ade=0; 


Zo 


das ist nichts anderes als das Integral von (gu), +0,= 0, also von der 
zweiten Gl. (15). 


$4. Eine nicht-lineare partielle Differentialgleichung zweiter Ordnung 
und ihre Integration durch „Aufspaltung“ 

Zweck der folgenden Betrachtungen ist, zu zeigen, daß die 
Gl. (17) eng mit einer nicht-linearen Differentialgieichung zweiter 
Ordnung zusammenhängen, deren Behandlung bei dem Zusammen- 
hang (2) soweit möglich ist, als die physikalische Fragestellung geht. 
Wir setzen 
(19) 
und erhalten aus (17): 


(20) 


m 
Daher setzen wir’): 
a=8); 
Wegen S, =D, folgt weiter: S=W,; U=W,. Ist nun: 
(21) 


so folgt: 
22 17 


1) Ein solches Eliminationsverfahren hat E. R. Neumann zuerst auf die 
hydrodynamischen Gleichungen angewendet (Crelles Journ. 132. S. 189). 
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als die gesuchte Die physikalischen Größen v, u 


sind aus W durch Differentiation zu finden: Aue Et 


(23) = =W a=W,, 

Es kommt also auf die Integration von (22) an; ist eine Lösung 
bekannt, so kann nach dem Verfahren der Gl. (23) und (18) das 
physikalische Verhalten in grundsätzlich einfacher Weise aus der 
Kenntnis von W abgeleitet werden. Die Integration (18c) ist nicht 
einmal nötig, denn aus (23) folgt: 7 = W,. Für den Zusammenhang (2) 
ergibt sich aus (22): ; 

a? 
24) Wie + 


n(W. 


m 


Wir wollen die Frage jetzt so stellen: Gegeben ist nur die 
Gl. (24); wie kann man sie integrieren? Ein Weg dazu ist in den i 
Arbeiten I und II vorgezeichnet; die jetzigen Ausführungen geben =| 
nichts grundsätzlich Neues dazu. Ich bringe sie, weil man der G1.(24) 
wohl nicht ansehen kann, auf welch sonderbare Weise sie integriert 
werden kann. Im Grunde ist das Folgende nichts anderes als ein 
Rückwärtsgehen des Weges von (24) zur Darbouxschen Differential- 7 ü 
gleichung. Man differentiiere (24) nach m; so wird: er 


Man führe z, = Z (es ist bis auf einen kon- 


n—1 


stanten Faktor das § der früheren Arbeiten) als neue Funk- 
tionen ein und bilde die Differentialquotienten von @ und = 
nach m und t: 


Daraus folgen unter Beachtung von (25) die Beziehungen: 
n+l 


1 1 
n+ n+ 


Das sind natürlich wieder die hydrodynamischen Ausgangsgleichungen, 
nur ist statt o die zweckmäßigere Veränderliche z eingeführt. Grund- 
sätzlich kann man jede durch Elimination entstandene Gleichung 
wie [(24)] natürlich wieder auf Systeme einfacherer Gleichungen 
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zurückführen. Aber es scheint mir lehrreich, den Weg von (24) aus 
einmal in der umgekehrten Richtung zu gehen. 

Die folgenden Schritte entsprechen denen in der Arbeit I; man 
betrachtet m und ¢ als Funktionen von uw und z und bekommt: 

n+1 n+1 
(27) m+Az"*-!t =0; m, + Az*-1t =0. 
n+1 

Jetzt setze man: m=V,; — Az"-1t=V, und erhält die Dar- 
bouxsche Differentialgleichung: 


V s+1 V, 


Damit ist die nicht-lineare G1.(24) oder (25) auf die lineare Gl. (28) 
zurückgeführt. Über die Integration von (28) ist in I und II ge- 
n+1 
sprochen; man bekommt schlieBlich m und — A z"~! tals Funktionen 
von z und u. Das bedeutet: (25) ist so integriert worden, daß m 
und ¢ als Funktionen von Z, und Z, bekannt sind. Beim Vergleich 
mit der Integration der Darbouxgleichung in II ist ein Unterschied 
in der Bezeichnung zu beachten: der konstante Faktor vor V,/z in (28) 
ist etwas anders aus n zusammengesetzt als in der Gl. (5), II. 
Ein sehr einfaches Beispiel zur Zerlegung einer Differential- 
gleichung in Gleichungen niedrigerer Ordnung ist die Aufspaltung 
der gewöhnlichen Wellengleichung (6a). Wir bekommen dabei auch 
eine Ableitung der d’Alembertschen Form der Lösung, die nicht 
umständlicher ist als die bekannte Ableitung, die von der Erfüllung 
der Potentialgleichung durch analytische Funktionen komplexen 
Arguments ausgeht. Zur Aufspaltung setze man U,= R, U,=8 
und erhält die zwei Gleichungen: 


(29) R,=8,; 8,= R,. 
Durch Addition und Subtraktion folgt daraus: ine 4 


(30) (R+ eS), = e(R+ eS); ®=1 


Diese Gleichung kann sofort integriert werden und gibt wegen der 
Bedeutung von R und S: 

(31) 

jede der Gl. (31) stellt ein erstes Integral von (6a) vor, ein auf den 
ersten Blick überraschendes Ergebnis; [aus der d’Alembertschen 
Lösung (6c) kann es natürlich leicht gewonnen werden]. Durch 


Addition und Subtraktion findet man U, und U, einzeln und durch 
Integration die Lösung (6c). 


4 gie 
Di 
an 
y 
da 
siit 
| 
qu 
wil 
unt 
T= 
sin 
Dif 
dal 
der 
a@al-, 
ähr 
Gl. 
gle 
oF 
in i 
x 
; das 
vie 
Ter 
— = 
laut 
(32: 
die 
| 
(321 
die 
anes 
(32: 


K. Beche rt. Dif beeiitiölgkeichen ngen n der we in Gasen 369 


In diesem Zusammenhang kam mir die Frage in die Quere, 
wieviel willkürliche Funktionen ein Integral von partiellen Differential- — 
gleichungen haben muß, damit es ein allgemeines Integral vorstellt. 
Die Frage läßt sich leicht an Hand der Taylorschen Reihe be- 
antworten. Man hat zu untersuchen, wieviel Glieder der Reihe Er 
beliebig vorgegeben werden können und welche durch die Differential- 
gleichungen und ihre Ableitungen bestimmt sind. Ein Beispiel wird 
das klarmachen; wir nehmen wieder die Wellengleichung (6a), grund- 
sätzlich ist die Überlegung auf beliebige Differentialgleichungen an- 
wendbar. 

Durch (6a) sind bestimmt: U,, und alle seine Difierential- 
quotienten nach gz und t; alle übrigen Glieder der Taylorreihe sind 
willkürlich wählbar. Entwickeln wir etwa in der Nähe der Stellez=0, 
t=0, so sind demnach willkürlich: U(0,0), U,(0,0), U,,(0,0), U,,,(0,0) 
und alle höheren Differentialquotienten von U nach z an der Stelle 
z=0,t= 0; daher ist die Funktion U (z,0) willkürlich; außerdem — 
sind beliebig wählbar: U,(0,0), U,,(0,0), U,,,(0,0) und alle höheren 
Differentialquotienten von U, nach x an der Stelle 2=0, t=0; 
daher ist die Funktion U,(z,0) willkürlich. Damit sind alle Glieder 
der Reihe erschöpft. Die allgemeine Lösung von (6a) muß daher 
zwei willkürliche Funktionen enthalten. 

Ich erwähne noch, daß man für Zylinder- und Kugelwellen ein 
ähnliches Eliminationsverfahren verwenden kann wie das zwischen 
Gl. (20) und (22) angewandte; es ergibt sich so eine Differential- — 
gleichung, die komplizierter ist als (22). a 


$5. Über das Problem der Ausbreitung ebener Wellen 
inidealen Gasen, wenn der Zusammenhang (2) nicht vorausgesetzt wird 


Wir schreiben die Beziehungen sogleich für den Fall hin, daß 
das Gas Reibung und Wärmeleitung zeigt. Volumkräfte sollen nicht _ 


wirken; die innere Energie pro Gramm bezeichnen wir mit U, Br mo 


Temperatur mit T, die Reibungskonstante mit u, die Wärmeleit- 
fähigkeit mit A. Die Kontinuitätsgleichung (Erhaltung der Masse) 
lautet: 


(32a) +ou,=0; 


die Bewegungsgleichung (Impulssatz): 


du 
(32b) a 


die Energiegleichung: 
d ut at = 
(820) ea (U+4 
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u und A sind dabei als konstant vorausgesetzt. Als Zustandsgleichung 
haben wir: 
(32d) p=oCT; 


und als Ausdruck fiir die innere Energie: 

C 
(32e) Um T + const; 


x ist das Verhältnis der spezifischen Wärmen Cy/Cy Auch hier ist 
es zweckmäßig, m,t statt z,t als unabhängige Veränderliche zu ver- 
wenden. Wie in den vorigen Paragraphen bedeuten die Querstriche, 
daß die Größen auf die Variablen m,t umgerechnet sind. Statt o 


a 1 
nehmen wir das spezifische Volumen v = — als Variable und er- 


Q 
halten an Stelle von (32a—c): er 
CT 
(33¢) T+ +),+ (c uülg?), — = 0. 


Ein Weg zur Untersuchung der Gleichungen ist wieder das Elimi- 
nationsverfahren; wir wollen es hier nicht im einzelnen durchführen, 
nur erwähnen, daß es auf eine nicht-lineare Differentialgleichung 
fünfter Ordnung führt!) und daß die Transformationsmethode von 
$ 2 auch hier geeignet erscheint, spezielle Lösungen zu verschaffen, 
Man kann die Transformationsmethode aber auch direkt auf die 
Gl. (33) anwenden. 

Im folgenden beschränken wir uns auf den Fall 4 = 0, u=(0; 
dann haben wir es wieder mit ebenen Wellen in Gasen ohne Reibung 
und Wärmeleitung zu tun; der Unterschied gegen die frühere Be- 
handlung in den Arbeiten I—III und die Überlegungen in den 
vorigen Paragraphen dieser Arbeit ist der, daß wir einen Zusammen- 
hang p = f(o) nicht voraussetzen. 

Durch Einsetzen von (33a,b) in (33c) folgt: 


T, 
T + (x—1) ; = 0, 
was sich integrieren läßt und gibt: 


(34) To! = f(m); 

f ist eine beliebige Funktion von m. Das bedeutet den allgemeinen 
adiabatischen Zusammenhang: 

(343) p = Cf(m)or, 


1) In & und ? geschrieben steht sie in der Dissertation von W. Kir 
hoff, Marburg 1931. 
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der sich von der sonst untersuchten Form (2) vor allem dadurch 
unterscheidet, daß an Stelle der Konstanten (in unserer früheren 
Schreibweise «a?/x) eine willkürliche Funktion von m steht. Physika- 
lisch bedeutet dies, daß die „Konstante“ p,/o, des Adiabatengesetzes 
von m abhängt, also für die verschiedenen Gasschichten verschiedene 
Werte haben kann. 
Das Eliminationsverfahren läuft nach dem Schema von Gl. tad 
bis (22) und liefert: 
(35a) a=W, 
und als Differentialgleichung für W: earere 
(35)) F (m) = 
F ist wieder eine willkürliche Funktion von m, die in einfacher 
Weise mit f(m) zusammenhängt. Ein Integral von (35b), das eine 
wilkürliche Funktion enthält, bekommt man durch die Unter- 
suchung der möglichen Ähnlichkeitstransformationen. (35b) erlaubt. 
die Transformation: 
W=aW; t=a ? t; m=m. 
Verlangen wir gemäß unserer in § 2 entwickelten Methode, daß 
speziell: 
„+1 


W’=W(m',t') = 3 ) = «W (m,t) 
sei, so folgt durch Differentiation nach a: 3° 


—1 
‘(x + 1) = 


also: 


(x + 1) 7 


W,. 


daher, wenn wir die Striche als jetzt unnötig weglassen: 
36) W=y(m)t 
Das ist in der Tat ein a von (35b), wenn “ee 2 
2(1— 
(«+ 
Da F willkürlich ist, ist auch g willkürlich und braucht nicht aus 
F berechnet zu werden. Man bekommt vielmehr mit Hilfe von 
(5a) und (36) unmittelbar Lösungen von (33) bei beliebigem g(m). 
Beachtet man noch, daß (33) außerdem die Transformationen t ->t + t,, 
">m+m, und u>u-+ 4 erlaubt — das letztere erkennt man 
durch Einsetzen von (33b) in (33c) — bei beliebigem konstanten t,, 
™,4, so wird unsere Lösung: 


F (m) = 


2 
| 
| 
. 
n 
4 
e 
): 
u 
“8. 
| 
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x—1 


9° ~ 
; galt +t) 4; 


2(x—1) 


‘= Im m (t + Im mit-+t,) 


2x 
2(x—-—1) 


Bei der Umrechnung auf die Veränderlichen x, ¢ kommt zu- 
nächst durch Integration von @ nach t: 


Beine beliebige Konstante. Weiter findet man, wenn h eine will- 
kürliche Funktion des Arguments rn j 


Enel 


2 mög) 


(38 a) At)(t+t,) **} 


t+ ty +4; t daß 
2x relat 


2(x—1), = 

= pe t+ 10) fydh; CT=pvr. weit; 

Die Wellen (38) verzerren ihre Form beim Weiterwandern. ofl 

Z.B. ist die Wanderungsgeschwindigkeit eines fest vorgegebenen bei d 
_u-Wertes bestimmt durch: 


die | 

dx (x + 1) 4(*—1) 
tatıv, 


also von u abhängig. Daher laufen die verschiedenen u-Werte ver- kon 
schieden schnell, die u-Welle verändert ihre Form. Keine Ver- uch 
zerrung würde es nur für «= — 1 geben, ein solcher x-Wert kommt weh 
aber bei Gasen nicht vor. Die Berechnung der Wanderungs- a 
geschwindigkeit der v-, p- und T-Wellen gibt verwickeltere Aus- pe 
drücke; auch diese Wellen verzerren sich. Dies Ergebnis war zu Sal 
erwarten, denn schon im Spezialfall F(m)=const, der in I und Il ‘inthe 
behandelt ist, hatten wir Wellenverzerrung (und damit die Möglich- fäche 
keit von Stoßwellen) erhalten. Einfü 

Ich freue mich, daß diese Arbeit in das Heft Aufnahme finden a 
wird, das dem 60. Geburtstag von Prof. Madelung gewidmet ist dark 
und bitte ihn, sie als nachbarlichen Glückwunsch anzusehen. werte 

Gießen, Institut f. theoret. Physik der Universität, 9. April 1941. und e 
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Von 8S. Flügge 


Die Eigenschwingungen eines Flüssigkeitstropfens werden für den all- 
gemeinsten Fall einer nicht-axialsymmetrischen Oberflächenschwingung be- 
rechnet. Die Lebensdauer eines solchen Zustandes gegen Abstrahlung von 
Energie wird nach der klassischen Elektrodynamik ermittelt is eine neue ae 
mögliche Auffassung der Kernisomerie darauf 4 

1. Einleitung‘ 

Seit einigen Jahren hat sich in zunehmendem Maße gezeigt, 
dad für viele Fragen der Kernphysik die Quantenmechanik eine 
relativ untergeordnete Rolle spielt und die klassische Mechanik 
weitgehend zu ihrer Behandlung ausreicht. Der Grund für dies 
unterschiedliche Verhalten zur Physik der Elektronenhülle liegt im 
wesentlichen darin, daß infolge der viel größeren Kernmasse der 
bei den Eigenschwingungen eines Atomkerns bewegten Materieteilchen 
die Nullpunktsenergie ziemlich klein ist im Verhältnis zur Schwin- 
gungsenergie der tiefsten Anregungsstufe. Es ist daher sicher quali- 
tativ, und mit einigem Recht auch quantitativ statthaft, den Atom- 
kern als ein Flüssigkeitströpfchen zu behandeln, dessen Energieinhalt 
sich zusammensetzt aus Volumenergie und Obertiichenenergie, 
wohinzu bei den schwersten Atomkernen noch in merkbarem Maße 
die elektrostatische AbstoBungsenergie der Protonen tritt. Die 
absolute Größe jedes Energieanteils kann dabei durch eine Material- 
konstante beschrieben werden, die für die Kernmaterie charakte- 
ristisch ist, z. B. die Größe der Oberflächenenergie durch die Ober- 
flichenspannung 7. Die Aufgabe der Kernphysik zerfällt durch die 
Einführung solcher Materialkonstanten, genau wie die Untersuchung 
makroskopischer Flüssigkeiten in zwei Teile: einen den Hilfsmitteln 
der klassischen Mechanik der Kontinua zugänglichen, der die Zahlen- 
verte der Materialkonstanten als Gegebenheiten a priori voraussetzt, 
und einen zweiten, tiefer dringenden Teil, der eben diese Konstanten- 
werte zu erklären strebt und dabei mit dem scharfen Werkzeug der 
Quantenmechanik an Fragen der Struktur herangeht. Wir wollen 
uns hier auf die erste Fragestellung beschränken, bewegen uns also 
= im Kreise der klassischen Mechanik und RN Es 
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. on meiner Ansicht nach zu den schönsten Ergebnissen der 
kernphysikalischen Entwicklung der letzten Jahre, daß sich immer 
klarer die Möglichkeit herausgeschält hat, solche rein klassische Be. 
trachtungen für viele Fragen anzustellen. In dieser Hinsicht wäre 
wohl vor allem auf den großen Erfolg nachdrücklich hinzuweisen, 
den Bohr und Wheeler gehabt haben bei der Behandlung des 
Spaltungsprozesses der schwersten Atomkerne nach dieser Methode! 

Das Ziel der vorliegenden Arbeit ist ein doppeltes: Einmal 
soll ein möglichst geschlossenes Bild der Methode gegeben werden, 
um künftige Arbeit auf diesem Gebiet zu erleichtern. Dazu ist es 
notwendig, manche Dinge, die sich in der älteren hydrodynamischen 
Literatur versteckt finden, in einer für unseren Problemkreis passenden 
Form darzustellen, auf die Gefahr hin, hier und da bekannte Dinge 
zu wiederholen. Zum andern soll eine neue Anwendung dieser Be- 
trachtungsweise, unter Hinzuziehung von Methoden der klassischen 
Elektrodynamik, auf den Problemkreis der Kernisomerie gezogen 


Wir stellen uns den Atomkern im folgenden als Tröpfchen 
einer reibungsfreien, inkompressiblen Fliissigkeit vor. Wie weit die 
Voraussetzung der Reibungslosigkeit erfüllt ist, läßt sich nicht leicht 
übersehen. Die Inkompressibilität ist nahezu erfüllt, wie schon 
ältere Betrachtungen von Bethe?) zeigen, solange man nicht zu 
allzu hohen Anregungsenergien übergeht. Wir beschränken uns ferner 
auf diejenigen Schwingungszustände des Tröpfchens, die keinen 
Drehimpuls haben. Wir treffen damit eine Auswahl unter den mög- 
lichen Zuständen, über die weiter unten noch zu sprechen sein wird. 
Endlich wollen wir noch die Coulombsche AbstoBung der Protonen 
vernachlässigen. Das bedeutet, daß unser Modell nicht ohne weiteres 
auf Vorgänge bei den schwersten Atomkernen anwendbar ist; gerade 
hier haben aber Bohr und Wheeler bereits gezeigt, in welcher 
Weise man vorzugehen hat, um der Coulombenergie Rechnung zu 
tragen. 

Die mathematische Formulierung unserer Voraussetzungen führt 
zu folgenden Aussagen: Die Schwingungen des Tröpfchens mögen 
beschrieben werden durch den Geschwindigkeitsvektor » jedes Volum- 
elements dr. Dann ist rot » = 0, d.h. es existiert ein Geschwindig- 
keitspotential: 

(1) | b= grad 


3 
1) N. Bohr u. J. A. Wheeler, Phys. Rev. 56. S. 426. 1939. 
2) H. A. Bethe, Rev. Mod. Physics 9. S. S6ff. 1937. 
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Wegen der Inkompressibilität ist außerdem divo = 0; also befolgt 
® die Differentialgleichung 


(2) AP=0. 


Bezeichnen wir mit r, +, gq die Kugelkoordinaten eines Punktes im 
Innern der Flüssigkeit, so lautet die allgemeinste bei r = 0 reguläre 
Lösung dieser Gleichung 


(3) ®= Pr(cos 


Ist (R, #,%) ein Punkt an der Oberfläche des Tropfens, so muß sich 
R ebenfalls in Form einer Entwicklung nach Kugelfunktionen dar- 
stellen lassen: 


(4) + Sia, „ Pi" (cos 9) 


Zwischen den Koalfisienten @,,, und #,,, der beiden Entwicklungen 
besteht ein einfacher Zusammenhang. Es muß nämlich die Radial- 
geschwindigkeit ÖR/Öt eines Oberflichenteilchens gleich sein der 
Ableitung (0 ®/Or), — z; mithin 

In dieser Formel ist bereits R durch R, ersetzt, d. h. die Voraus- 
setzung eingeführt, daß die Schwingungsamplituden &, ,,< R, sind. 
Das bedeutet die Harmonizität der Eigenschwingungen. 


3. Berechnung der Eigenfrequenzen des Tropfens 


Um die Eigenschwingungen näher zu untersuchen, stehen wie 
bei allen mechanischen Problemen zwei Wege offen: Entweder be- 
ginnt man mit der Betrachtung des Kräftespiels, oder man geht 
unmittelbar vom Energiesatz aus. Während der erste Weg um seiner 
Anschaulichkeit willen von der klassischen Hydrodynamik bevorzugt 
wird, entspricht der andere mehr den Denkgewohnheiten des Atom- 
physikers und soll deshalb hier eingeschlagen werden’). 

Ist o die Massendichte der Kernmaterie, so ist die kinetische 
Energie der Schwingung 


xe f |grad dr, 


wobei die Integration über das ganze Volumen des Trépfchens zu 
erstrecken ist. Die potentielle Energie bei einer Vergrößerung der 


1) Der Kräfteansatz ist durchgeführt im Handbuchartikel von Gyemant, 
Handb. d. Physik VII, S. 366f. Unsere Rechnung geht inhaltlich darüber hinaus, 
da dort nur das unvollständige Problem rotationssymmetrischer Schwingungen 

"(m = 0) behandelt ist. 
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Oberfläche um den Betrag 4F wird y- AF. It C=R— Rk, die 
radiale Verschiebung eines Oberflächenpunktes und d2 = sin itd ddg 
ein Raumwinkelelement, so nimmt die potentielle Energie die Gestalt an: 


Setzen wir in Gl. (6) das Geschwindigkeitspotential ® aus Gl. (3) ein, 
so erhalten wir die Ausdrücke: F 


= } rsind d@ 
> Br + Deo (1 sin. P™+cos 
l,m \ dd sins 

cos? O@D 
(sin or + r rsin® 


—y dP," m 
PP» 


r 


sin 3S) 


m 
= D Bunt (I cos # P=— sin # = =). 


Bildet man aus diesen Ausdriicken und den komplex konjugierten 


1 (8% 4580) 


oy 


(ag 
+ 


6 
so sieht man unmittelbar ein, daß in den Vierfachsummen über 
l,l’, m,m’ nur die Glieder mit m = m’ zur kinetischen Energie bei- 
tragen. Nach Ausmultiplizieren erhält man den Ausdruck: 


1 I+V ( 


om m 2 \ 
‚pm m aP, dP; m 


In ganz ähnlicher Weise, nur einfacher, ergibt sich aus Gl. (4 
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/ Winkel ountitilions, bedienen wir uns der Hilfsformeln ): 


ap 

um die Ableitungen umzuformen. Die dabei iibrigbleibenden Glieder Pte 
mit cot + formen wir dann um mit Hilfe der aus G1.(10) folgenden 
Beziehung 


IM mct#P’ = m) (lL—m + 1) 4 pret, 


Es läßt sich dann unmittelbar die Orthogonalitäts- und Normierungs- ar 


(12) Q Pp” (9 


mit der sich nach längerer, aber einfacher Rechnung die Ausdrücke a 
ergeben: 


4a (+ m)! 
(13) Exin = 5 > 7 Bim 1(2141)- Foi? 


[ = — 9 é 
(14) Epot 57 > >| wer d—1)(+2). 
m 

Unter Heranziehung der Gl. (5) zwischen den Koeffizienten «, „ 
und £,,, entsteht so der Ausdruck fiir die Energie: 

r 


(+ m)! 


mit der Abkiirzun 


Aus dem Energieausdruck (15) folgt in bekannter Weise zunächst 
die Unabhängigkeit der einzelnen, durch Wertepaare !, m charakte- — rn 
risierten Eigenschwingunger, sowie ihr sinoidaler Charakter in der a 

Zeit mit der von m unabhängigen Frequenz w,. Wir dürfen also = 
setzen 


(17) = Ay, COS 4 


—_ 


1) Für Hilfsformeln über Kugelfunktionen wurden benutzt: E. Madelung, © 
Die math. Hilfsmittel des Physikers, 3. Aufl.; Jahnke-Emde, Funktionen- fees 
tafeln, 2. Aufl. sowie H. Bethe, Hab. d. Physik 2. Aufl. Bd. XXIV. 1. S. 551ff. a oe 
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die Energie der Schwingung wird dann 


(l + m)! 
+ 2) a; 


. 
l,m 


4. Einige Anwendungen auf den Atomkern 


Wir wollen die Ergebnisse des vorigen Abschnitts dazu be- 
nutzen, einige Betrachtungen über die angeregten Zustände der 
Atomkerne anzustellen. Aus dem Tröpfchenmodell folgt zunächst 
qualitativ die Möglichkeit, die Terme eines Kerns nach steigender 
Energie in drei Gruppen einzuteilen, ähnlich wie man bei den Eigen- 
schwingungen eines Moleküls Rotations-, Schwingungs- und Elek- 
tronenterme unterscheidet. Diese drei Gruppen sind: ; a 


a) Rotationsterme 


Es genügt infolge des großen Trägheitsmoments J bereits eine 
geringe Energiezufuhr, um den Kern in Rotation zu versetzen. Wegen 
der großen Langsamkeit einer solchen Rotation dürfte der Kern 
wie ein starrer Körper rotieren, d. h. die möglichen Drehimpulse 
sind Jh und die zugehörigen Energiestufen 


2 
2) 
Hierin kann man J = 5 MR,? einsetzen; mit der Formel 


(20) R, = 1,4. 107°. cm 

5 
erhält man dann E,= 274 J(J-+1)MeV; für J - A =100 
führt das zu einer Anregungsenergie von 25 keV. 


b) Oberflichenschwingungen 
Bei etwas größerer Energiezufuhr ist der Kern auch in der 
Lage, Deformationsschwingungen auszuführen mit den im vorigen 
Abschnitt berechneten Eigenfrequenzen @,, d.h. also mit den An- 
regungsenergien E,=4w,. Nach Gl. (16) lautet die Termformel 
hierfür 


Setzt man den empirischen Wert 427 R,? = 14- 4": MeV für die 
Oberflichenenergie ein, so erhält man numerisch 


MeV. 
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Der Zustand mit /=0 ist der Grundzustand; eine Schwingung 
mit {= 1 bedeutet lediglich eine starre Parallelverschiebung des 
ganzen Kerns. Sie ist also physikalisch sinnlos, was sich darin 
äußert, dab ©, = 0 wird. Die Schwingungen mit != 2 sind die 
tiefsten echten Schwingungszustände, bei denen eine ellipsoid-ähn- 
liche Verformung eintritt. Die Anregungsenergie dieses 5-fach ent- 
arteten Zustandes beträgt bei A = 100 etwa 2,8 MeV. Das ist 
viel größer als E,. Es wird also eine ziemlich weitgehende Ent- 
kopplung der Rotations- und Deformationsschwingungen eintreten. 
Übrigens verfügen wir kaum über experimentelle Möglichkeiten, 
Terme mit hohen Rotationsquantenzahlen anzuregen. 


ec) Dilatationsschwingungen 


Die Anregungsenergie der Schwingungen, bei denen die Kern- 
dichte nicht konstant bleibt, wird um so größer sein, je geringer 
die Kompressibilität des Kerns ist. Diese abzuschätzen, ist schwierig; 
sie zu berechnen, setzt eine nähere Kenntnis von Strukturfragen 
voraus, die wir heute noch nicht annähernd besitzen. Bethe!) hat 
versucht, die Kompressibilität auf Grund des Fermi-Gas-Modells 
abzuschätzen; er findet, daß die Dilatationsschwingungen erst bei 
ziemlich hohen Anregungsenergien (= 10 MeV) spürbar werden. Hier 
ist es aber wohl nicht so, daß man eine weitgehende Entkopplung 
von Deformations- und Dilatationsschwingung hat: Dieser Punkt 
trägt daher viel bei zur Komplikation der theoretischen Deutung 
der Kernspektren. — 


Die Berechtigung unseres Modells können wir noch in anderer 
Hinsicht überprüfen. Wir wollen nämlich nachsehen, ob die Ampli- 
tude a, der Schwingungen wirklich klein ist gegen R,. Wir dürfen 
uns dabei auf den ungiinstigsten Fall m =0 beschränken. Setzen 
wir in Gl. (18) für E,, den Wert von Aw, nach (21) ein, so er- 
halten wir 
(22) a; ) V: A? 


\ 


Setzt man hier Zahlenwerte ein, so erhält man ne 


a \* _ 1) 
(x) “Va — 1042) 


und die Größenordnung R= 0,16 für A = 100. 


1) H. A. Bethe, Rev. Mod. Physics 9. S. S6ff. 1937. 
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5. Strahlungstheorie für das Tröpfehenmodell des Atomkerns 


Da ein Atomkern die Ladung Ze trägt, muß er bei einer 
Schwingung nach der klassischen Elektrodynamik Energie abstrahlen. 
Diese Strahlungsdämpfung ist das korrespondenzmäßige Äquivalent 
der Übergangswahrscheinlichkeit aus dem angeregten Zustande der 
Energie Aw, in den schwingungslosen Grundzustand. Es ist unser 
Ziel, auf diesem klassischen Wege die mittlere Lebensdauer einer 
angeregten Deformationschwingung w, eines Atomkerns zu berechnen, 

Es sei o die Ladungsdichte, die im Innern des Atomkerns 
konstant angenommen werden darf. Im Volumelement dr an der 
Stelle r (Kugelkoordinaten r, :+, g) befindet sich dann die Ladung o dr. 
Wir betrachten einen Aufpunkt r, sehr weit außerhalb des Tropfens; 
seine Kugelkoordinaten seien r,,&,,%,. Netzen wir noch 1,—t=3, 
so werden von einer Bewegung der Ladung odr die Feldstärken 
an der Stelle en erzeugt: 


d& odı 3 5 
em a „ex 
oat 
| dQ = 
Dabei bedeutet ¢ nun nicht BEN die zweite Ableitung nach der 


Zeit, sondern es ist noch die Retardierung zu KENN 
Wollen wir z. B. die Feldstärke im Aufpunkt zur Zeit != + 


berechnen, so müssen wir in ¢ die Zeit ?— einsetzen. Ist die 
Bewegung insbesondere eine harmonische Schwingung, so müssen 
wir also in Gl. (23) einführen: 

r,—'8 
t = fonit. + acoso (t+ = 


| und 


| =— aw? w(t =). 


Nennen wir den von den V arena r, und r eingeschlossenen W inkel 9, 


so wird 
— ro — — ‘) . 


Macht man = < = durch geeignete Wahl des Aufpunktes (r, > & 


beliebig klein, so genügt das erste Entwicklungsglied: 


(25) = == . 2080 = cos 0. 

Wie man sieht, wird dies erste Glied unabhängig vom Aufpunkt. 
Es gibt die Phasendifferenz der von verschiedenen Punkten des 
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Kernes ausgegangenen Wellenzüge am Aufpunkt an. Da“ = = = 5 


mit der Wellenlänge zusammenhängt, schreiben wir ews 


2 


(26) — 8) = cos 


Diese Größe können wir nicht beliebig klein machen; bei Problemen 
der Kernphysik ist sie aber fast immer klein genug, um Potenzreihen- 
entwicklungen zu gestatten, mit deren erstem von Null verschiedenem 
Gliede man sich begnügen kann. Wir schreiben deshalb 


cos @ (1 + - =) = cos (ot + cos 0) 


fe} 
(5) cos Ox...| 


2 
— sin E -cosO — - (5) c08? O 


+ OF... 


Setzen wir an t nach Gl. (24) und (27) in Gl. (23) ein, ersetzen 
noch, was ja erlaubt ist, den Vektor $ durch rt, und integrieren 
über alle Volumelemente dr innerhalb des Tropfens, so entstehen 
unter Einführung der Abkürzung 


die Feldstärken 


2 
cr, 


im Aufpunkt. Um den gesamten Energiestrom zu erhalten, der den 
Flüssigkeitstropfen in Form von Strahlung in der Zeiteinheit ver- 
läßt, bilden wir den Poyntingvektor!) 


(30) =~ Exh) = 


der radial nach außen weist, wie es sein muß. Insgesamt wird je 
Zeiteinheit vom Tropfen die Energie 


(31) s-[f 2, = - (x, 


Anchr,! 


2972 2 
— (r, M4 r,, 


5 t 
4nc*ry 


abgestrahlt. 


1) Die den vektoriellen Umformungen zugrundeliegenden Hilfsformeln 
finden sich bei E. Madelung, a. a. O. 
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Unsere erste Aufgabe wird es sein, den in Gl. (28) definierten 
Vektor X zu berechnen, den wir nach Gl. (24) und (27) ausführlich 
schreiben können: 
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| - [dr a(r, 


{cos w t- 1-3 


(32) ; 
-\5 
Für den Amplitudenvektor a greifen wir auf die Gl. (1) und (3) zurück, 
wobei wir uns jetzt auf eın Glied der Summe in Gl. (3) beschränken 
können, da wir bereits wissen, daß die einzelnen Glieder unabhängige 
Eigenschwingungen entsprechen. ‘Wir haben also (abgesehen vo 
der Retardierungsphase) einerseits 


®, = P,r!P,"(cos = — o,b,r! P™(cos sin w, t, 
wodurch wir die Größe b, definieren wollen, anderseits nach Gl. (24) 
b=grad®, = — w,-asinw,t. 

Der Vergleich beider Ausdrücke führt auf i 
(33) a = b,grad (r! Pe'™?), 


Die Größen b, können mit Hilfe von Gl. (5) und (17) durch die 
a, ausgedrückt werden: 


(34) a, = b,-1Rj-?. 


Aus den Gl. (31) bis (34) haben wir jetzt alle weiteren Schlüsse’ zu 
ziehen. 
Der Vektor a hat die Komponenten 


Aug 


dP,” 
(isin + + cos = 


a, — va, =b = (sin + + cos + — 


a,=b = byr'-1eime (Leos — sin )- 
Diese komplizierten Ausdrücke lassen sich auf je eine Kugelfunktion 
der Ordnung (/— 1) reduzieren. Hierzu muß man zunächst unter 
Verwendung der Hilfsgleichungen (10) die Ableitungen wegschaffen 
und alsdann die Formeln anwenden 
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(21+ P* 


m m 1 
21+1)sin#P”=P — p**', 


| (27+ 1)sin P™ = — m+ 1)(l— m+ 


+(+mP,_, 


(35) 
l+1 


d+ mi+m— 


Es ergeben sich so die Gleichungen 


= Ob, m) +m— (iF), 
= br” m)e'"? P®_,(#). 


Die Berechnung der Komponenten von U läuft also darauf hinaus, 
Integrale zu berechnen, die jeweils über das Produkt einer Kugel- 
funktion der Ordnung (/— 1) mit einer Potenz von cos® zu er- 
strecken sind. Um das durchzuführen, ersetzen wir die Potenzen 
von cos@ durch Legendresche Polynome; es gelten nämlich die 
Identitäten: 


10 
P — P,(9)+ 37 P,(0)+ = 


Führen wir diese Ausdrücke in Gl. (32) ein und ordnen nach 
Kugelfunktionen, so nimmt die geschweifte Klammer die Form an: 


P,(0)- 
| — sin wt P,@)| — 5 


+ coswt+ P,(9)- 


1 


— sin ot P,(0)| 445 


+cosot.P,(0)- 


; 


= 
F 
a | + | | 
Ä 
(36) a. 2 
| 
2 
N 
- 
. > | 
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i Wegen > <1 können wir uns in jeder der eckigen Klammern auf , 
U 
das erste "iea beschränken; es bleibt dann die Summe n 
- n+Dn r \n 1 (2 
(834) (ot (4) "73.5... Qn — 1) 
n=U 
i 
übrig, die wir als zweiten Faktor unter das Integral in GI. (32) zu / 
setzen haben. Die Anwendang der Normierungs- und Orthogonalitäts- d 
beziehung ') 
reduziert die Summe dann jeweils auf ein einziges Glied. Ersetzen 
wir noch b, gemäß Gl. (34) durch a, und führen an Stelle der 
Ladungsdichte o entsprechend der Beziehung 420 R,* = 3Ze die 
— Kernladung ein, so entstehen die Gleichungen 
(39) 1%, = A’ =k(l + "7 
U, = k(l + mie Pe 1 I 
in denen die Abkürzung 
3. .. In\ ( R,\!-1 
(40) k = 3Ze a, Ql +1)! sin (wt + 7) 
für den allen Komponenten gemeinsamen, winkelunabhängigen Faktor 
benutzt ist. 
Die Berechnung von S nach Gl. (31) erfordert die Kenntnis 
von |X|? und (r,W. Während man mühelos 
2 2 1 +2 2 all m+1\2 
‘ ; 1 2 2 m—1\2 2 m \2 
7 is 
bildet, sind bei a 


sin d, A + At) +17, cos J, AU, 


1) Vgl. E. Madelung, a. a. O., 8. 66, Gl. (30). 


Win 


r 


pts ps einige Umformungen auf Grund der Gl. (35) erforderlich, um zu dem 
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zu gelangen. In Gl. (31) treten bei Verwendung der Beziehungen (41) 
und (42) nur mehr Integrale vom Typus der Gl. (12) auf. Die Aus- 
rechnung ergibt in einfacher Weise 

Bilden wir das zeitliche Mittel S, ersetzen nach Gl. (21) und (22) die 
Amplituden a, unter Biiminisrang der Oberflächenspannung y, und 
dividieren endlich durch Ao,, so erhalten wir die Zahl der je Zeit- 
einheit emittierten Lichtquanten 


La 


ee was dasselbe ist, die Wahrscheinlichkeit für den TR des 
angeregten Zustandes in den Grundzustand in der Zeiteinheit: 


ho z\3¢-) 
(44) w =f,n'2 Me 


mit 
(45) 


Die Halbwertszeit des angeregten Zustandes hängt hiermit zusammen 


nach der Relation shale 
In2 


(46) | 


~~ “= 6. Anwendung auf die Kernisomerie 


Das erhaltene Schlußergebnis bedeutet wegen der Kleinheit 
des Verhältnisses von Kernradius R, und Wellenlänge 4 der emittierten 
y-Strahlung, daß die Lebensdauer eines mit der Quantenzahl / 
schwingenden angeregten Kernzustandes um so länger wird, je größer 
list. Die Gl. (44) zeigt dabei in ihrem ganzen Aufbau eine auBer- 
ordentliche Ähnlichkeit mit denjenigen, welche v. Weizsäcker!) und 
Bethe?) bereits früher zur Erklärung der Erscheinung der Kern- 
isomerie aufgestellt hatten. Während bei diesen früheren Formeln / 
aber stets die Bedeutung eines Drehimpulses hatte, so daß die 
Forderung nach hohen Drehimpulsdifferenzen für die tiefsten Terme 
aufgestellt werden mußte, handelt es sich hier um drehimpulslose 
Zustände, bei denen die Quantenzahlen / und m lediglich die An- 
zahl der Knotenlinien messen, mit denen die Kernoberfläche schwingt. 


1) C.F.v. Weizsäcker, Naturwiss. 24. S. 813. 1936. 
2) H. A. Bethe, Rev. of Mod. Phys. 9. S. 220ff. 1937, insbesondere 
Gl. (732) auf §. 226. 
ar der Physik. 5. Folge. 30. 
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Die hier gebotene Theorie geht an Geschlossenheit über die 
älteren Ansätze hinaus insofern, als es sich dort stets nur um 
einigermaßen rohe Größenordnungsabschätzungen handelt. Das er- 
kennt man schon daran sofort, daß die Berechnung einer Lebens- 
dauer bei mittelschweren Kernen nach den Formeln v. Wei zsiickers 
und Bethes für !=4 oder 5 um einen Faktor 10% bis 105 ver. 
schiedene Ergebnisse liefert. Einen Versuch, diesem Übelstande 
abzuhelfen, stellt die Arbeit von Koyenuma!) dar, die aber ihrem 
Wesen nach nur einen ersten Anfang bedeuten kann. 

Unsere Formel ist streng in dem Maße, in dem korrespondenz- 
mäßige Bilder für den Atomkern erlaubt sind. Es hat daher einen 
guten Sinn, auch die Größe des Normierungsfaktors f, „ anzugeben 
und daraus Schlüsse zu ziehen. 

Man sieht aus Gl. (45), daß für m=0, also für rotations- 
symmetrische Schwingungszustände, die Lebensdauer am größten 
wird. Die sich dann ergebenden Zahlenwerte f,, sind im folgen- 
den für einige | angegeben: 
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0,250 

1,09 - 
1,68 + 10-5 
1,66 - 10-7 
1,14 - 10-® 
5,81 - 10-18 


l 
2 
3 
4 
5 
6 
7 


Schon dieser Normierungsfaktor ergibt also durch seine bemerkens- 
werte Kleinheit eine mit steigenden / rasch anwachsende Lebensdauer. 
Die Gl. (44) enthält die gleichen Faktoren wie die Formeln 
v. Weizsäckers und Bethes: e?/kc, die mit 21 steigenden Potenzen 
von R,/# und den Faktor w, der dem ganzen Ausdruck die 
Dimension einer reziproken Zeitgibt. Der bei Bethe verschwundene 
Faktor Z? der ursprünglichen Weizsäckerschen Formel ist wieder 
vorhanden. Außerdem ist aber noch der dimensionslose Quotient 
der Anregungsenergie Aw, zur Ruhenergie Mc? des Atomkerns auf- 
getaucht, der davon herrührt, daß in die Deformationsschwingungen 
die Trägbeit der Kernmasse eingeht. Dieser Quotient ist natür- 
lich sehr klein, so daß er — ähnlich wie f,„ — abermals eine 
Herabsetzung des erforderlichen !-Wertes bedeutet. 
Im folgenden sind für einen mittelschweren Kern (R, = cm, 
Z=50, A=120) aus dem Gebiet, in dem nach den experimentellen 
1) N. Koyenuma, Ztschr. f. Phys. (im Erscheinen). a ee 
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Befunden die Kernisomerie am häufigsten auftritt, die Halbwerts- 


zeiten angeregter Deformationsschwingungen fir m=0 und !=3, 4 
und 5 bei verschiedenen Anregungsenergien E=4w, numerisch 
angegeben: 


Halbwertszeit 7 für 
l=4 


280 Jahre 
68 d 

6,3 h 

1,5 min 
0,06 sec 


Es scheint hiernach, als sollte man diese Betrachtungsweise be- 
sonders zur Erklärung derjenigen Fälle von Isomerie heranziehen, 
bei denen die beobachtete Anregungsenergie verhältnismäßig groß 
it. Doch bedürfen alle ins einzelne gehenden Fragen noch näherer 
Untersuchung, weshalb hier auch nicht mehr darauf eingegangen 
werden soll. 


Berlin-Dahlem, Kaiser - Wilhelm - Institut für Chemie, im 


April 1941. 


gegangen 10. April 1941) 


TER 
ch 


Br 
= 
| 0,1 0,003 „ 8000 Jahre 
| 
0,2 8 ” 
| 0,5 7,2 h 
| 1 25 sec 
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Reflexion elektromagnetischer Wellen 
an einer inhomogenen Schicht a 
nach der Wentzel-Kramers-Brillouin- Methode = 


(Mit 4 Abbildungen) 


Einleitung und Zusammenfassung 


Als Grundlage von Experimenten, in denen die Reflexion oder 
der Durchgang von elektromagnetischen Wellen durch inhomogene 
Schichten eine Rolle spielt, ist es wünschenswert, eine Formel für 
das Reflexionsvermögen zu besitzen, die möglichst weitgehend vom 
speziellen Verlauf des Brechungsindex in der inhomogenen Schicht 
_ unabhängig ist. Um eine solche Formel abzuleiten, benützen wir 
E im folgenden die Wentzel-Kramers-Brillouin (W.K.B.)-Methode°), die 
von ihren Anwendungen in analogen wellenmechanischen Problemen 
a bekannt ist. Volle Analogie mit den dort behandelten Problemen 
besteht in der Optik allerdings erst für die Totalreflexion, während 
wir uns hier mit senkrechtem Einfall begnügen. Unser Ziel ist 
damit enger gesteckt als das von R. Gans’), der bei seiner viel 
umfassenderen Behandlung schon im wesentlichen die gleiche mathe- 
_ matische Methode angewandt hat. Dagegen haben wir höhere 
Näherungen und die Interferenz zweier Sprungstellen des Brechungs- 
index oder seiner Ableitungen in Betracht gezogen. 

Ein Mangel, der allen bisherigen Behandlungen einschließlich 
der unsrigen anhaftet, besteht darin, daß nach ihnen das Reflexions- 
vermögen eines inhomogenen Mediums, dessen Brechungsindex selbst 
stetig und in allen Differentialquotienten stetig (z. B. vom Werte 
nı bei z =— oo nach dem Werte nq bei z=-++ oo geht), nicht be- 
rechnet werden kann. Reflexionen finden danach nämlich nur an 
Sprungstellen des Brechungsindex oder einer seiner Ableitungen statt, 
Nach einer kurzen Diskussion der Wellengleichungen für € 
und § in einem inhomogenen Medium (§ 1) stellen wir in § 2 für 

das Reflexionsvermögen einer inhomogenen Schicht eine exakte 
Formel auf, die wir durch Einführen der W.K.B.-Methode in eine 
‘ Näherungsformel verwandeln, welche nur noch den Brechungsindex 


und seine Ableitungen enthält. Sie ist gültig, solange (4 | 
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ist und bei entsprechenden Voraussetzungen über die höheren Ab- 
leitungen des Brechungsindex. In $ 3 vergleichen wir diese Naherung 
mit einem von M. Schlick!) exakt gelösten Spezialbeispiel, in 
welchem der Brechungsindex in der inhomogenen Schicht den Verlauf 
n(2) = 2 besitzt. In § 4 findet man das Reflexionsvermögen von 


Schichten, an deren Grenzen entweder der Brechungsindex selbst, 
oder seine erste, oder seine höheren Ableitungen eine Sprungstelle 
besitzen und in $5 wird die Frage besprochen, wann das Reflexions- 
vermögen einer inhomogenen Schicht das gleiche funktionelle Aus- 
sehen besitzt wie das einer homogenen Schicht eines mittleren 
Brechungsindex. 

Eine ausführlichere Diskussion und Anwendungen der er- 
haltenen Formeln werden in der Dissertation des einen von uns 
E. Menzer) erscheinen. 


§ 1. Die Wellengleichungen in inhomogenen Medien 


Wir betrachten den Durchgang von elektromagnetischen Wellen 
durch ein inhomogenes, nicht absorbierendes, isotropes Medium der 
magnetischen Permeabilität 1. Da in einem Medium mit veränder- 
licher Dielektrizitätskonstante 


divD = div(e ©) = ediv€ + (E, grade) = 0 


ist, treten in den Wellengleichungen fiir © und  Zusatzglieder 
auf, es gilt bekanntlich 


— 4€ — grad, grad In e) = 
) — 49 + [rot §, grad Ine] = 0. 


Solange jedoch © senkrecht grade steht, fällt das Zusatzglied in 
der Wellengleichung für E weg, solange € parallel grade liegt, ver- 
schwindet das Zusatzglied in der Wellengleichung für 9. 

Nehmen wir an, das Medium sei in der z-Richtung von z = z, 
ab inhomogen, ¢ = ¢(z) = n?(z). Im Falle senkrechter Inzidenz einer 
ebenen Welle z. B. E, = H, behält die Schwin- 
gungsgleichung für E, auch im inhomogenen Medium also ihre ge- 
wohnte Form 
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Um das von der Inhomogenität von & in der Wellengleichung für § 
herrührende Glied braucht man sich nicht zu kümmern, denn aus 
einer Maxwellgleichung folgt nach Auffinden der Lösung für E, 


wobei sich 
(4) 
ergibt. 


Auch bei schiefem Einfall fällt das Zusatzglied in (1) für die Komponente 
von €, welche senkrecht zur Einfallsebene liegt weg. Für diese Komponente 
führt also dasselbe Verfahren zur Lösung. Für die parallel zur Einfallsebene 
schwingende Komponente von € löst man am einfachsten Gl. (2) für die dann 
senkrecht zur Einfallsebene schwingende Komponente von § und erhält & aus 
der anderen Maxwellgleichung. Gl. (2) erhält dann zwar noch ein Zusatzglied, 
aber $/n befriedigt wieder eine Gl. (3) analoge Schwingungsgleichung, in der 
kein Glied mit erster Ableitung auftritt und in welcher an Stelle von n?o!/e 


n? 1/d 
der Ausdruck a + ; (7 In tritt. 


Br § 2. Das Reflexionsvermögen bei senkrechtem Einfall 


Wir bestimmen das Reflexionsvermögen in folgendem Fall: Aus 
einem von z =—oo bis z = z, homogenen Medium I des Brechungs- 
index ny fällt senkrecht eine zwischen z= 2, und in- 
homogene Schicht III des variablen Brechungsindex n (2) eine elektro- 
magnetische Welle. Von z=2, ab sei das Medium II wieder 
homogen vom Brechungsindex ny. Zur Lösung der Schwingungs- 
gleichung (3) bzw. 


3) 


in der inhomogenen Schicht III benützen wir die W.K.B.- Näherungs- 
methode. Für die beiden Lösungen von (3’) setzen wir danach 


2 z 
- | f(@)dz [mr 
F,@)=e * F_()=e 


worin f(z) = p(z) + iq(z) eine komplexe Funktion (p, q reell) und 2, 
einen beliebigen Punkt in III bedeuten. f* ist das Konjugiert- 
Komplexe von f. Beim Einsetzen der ersten Rn. (5) in (3’) entsteht 
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Die Gl. (3”) ieee wir durch re von f nach wachsenden 


Potenzen der Vakuumwellenläng 
wodurch Gl. (3”) in 
| Jo" —n)+ Qf, — ify) + + fi? — tf’) 
| + (2f,f,+ =90. 


Wir erhalten die sukzessiven FERNER aus (3°) durch Nullsetzen 
der Koeffizienten fortschreitend höherer Potenzen von c/w 


3” 


13 (n”)? | 99 n’?n” 
2 m: 2 


Allgemein gilt fiir die nte Näherung f, die Rekursionsformel 
a 2 fofn + 2hi fn -1 


fi für gerade n 


n/2 


| +t 


(8) 
+2fafn-2+ 


für ungeraden. 


Das negative Vorzeichen der fos 
zur zweiten Lösung F_(z) in (5). 


Realteil p und Imaginärteil g der Funktion f(z) hängen in einfacher 

Weise miteinander zusammen. Es ist bekanntlich 7 
BEA wre 

Dr: 
wie man durch Einführen von p und gq in Gl.(3”) erkennt. Da p nur die 
geraden fa, fa; q nur die ungeraden Näherungen fi, fa, fs; --- ent- 
hält, lassen sich die geraden aus den ungeraden berechnen. Für die prak- 
tische Auswertung der Rekursionsformel (8) entspringen daraus jedoch keine 
Vorteile. Aber der Lösungsansatz (5) läßt sich unter Benutzung dieser Tat- 
sache mit einer einzigen reellen Funktion p(z) in der bekannten Form 
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schreiben. Der Lösungsansatz (5) versagt, wenn der Realteil p(z) der Funk- 
tion f(z) verschwindet, denn dann stellen F, und F_ nur eine und dieselbe 


Lösung dar. Einsetzen von p und q in Gl. (3”) ergibt t 
2 2 2 + 14 + 1 9 c 0 
(2220 - <p = 0. 


Wenn sowohl p wie q existieren, so erhält man durch Trennen von Real- und 
Imaginärteil 

9 2pq- ee oder Gl. (9). 


-_ 


Verschwindet aber p und mit ihm p’, da q endlich bleiben soll, so entsteht 
für q folgende Differentialgleichung 


24? — ‘=(@. 
(9 


In den homogenen Schichten I und II besitzen die Lösungen die 
Form von ebenen Wellen mit konstanten Amplituden: Im Medium I 
sei die einfallende bzw. reflektierte Welie 


8,2 sot 
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| Ein A 
(11) 
| E' = A, 


im Medium II sei (die hier einzige) durchgelassene Welle z 
(12) E°=A,e =—n,,4,e : 
/ y 2 23 


z 


Im inhomogenen Medium III benutzen wir die Lösungen (5) für die 
einfallende bzw. reflektierte Welle 


Zur int C iat 

EM = A,F A, F, pe, 
H, =—1t—A,F_ (ge . 


Nach Berücksichtigung der Grenzbedingungen beim Übergang vom 
h Medium I nach III: 
y 


> < > < > < > < 
+E) =E™"+E™ und +H", 


und beim Übergang vom Medium III nach II: 
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erhält man für das Reflexionsvermögen __ 


2 2 
. € > e 
ny Nyy W, +i— (ny Wyo _ N] W,)+ WwW, 


| Ny Nyy Wy Wyo 


= F, (@,)F_(z,) —F_(z,)F, (2), 
a= F_(@)F, (2) — F, (2,) F_(z,), 


= 2) F_ (z,). 


Für die Ausdrücke (15) ergibt sich unter Anwendung des Ansatzes (5) 
und dessen Ableitungen 


FL@= L e (p + iq F, 


2io 
(p, ecoss+q, Sin 8), 


— 


P (20) 
Vii Ps 

bedeuten. Mit den Werten (15) nimmt der Ausdruck (14) für das 
Reflexionsvermögen folgende Form an 


(14) _ 9 8in’s +g, cos*s+g,sin2s _ G2 +9, + 29, tg 8 
9, sin?8 + 9, cos* s +g, sin 28 9s + 9,tg?s + 29, tgs” 


ir 


Aye 
T 
| 
( 
. fa 
| 4 wart 
e 
| 
| 
bzw. 
und der Anw endung der Abkürzungen p,=p(2), p,= Pl), 
| 
(15) 
W, Pi Pah J) 81N + (P, 
: 
wobei s= f 1 ie in o sich wegkürzende Größe 
ie 
U 
J 
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wobei unter 
2 


Nr n, 9; — 2m, My, P, Py + + (P,P. +4, 9.) 
P; pos 


9. = (", "yy P,) + — = 9,—4N Po, 
| I; = mi P24. + (Py Pot 9,9) (Pi Pod) 


verstanden werden soll. Der Ausdruck (14’) fiir das Reflexions- 
co _ vermögen gilt exakt, solange p und q existieren, denn die Näherung (6) 
ist in ihm noch nicht angewandt worden. 

Unter Annahme, daß der Nenner in der ersten Form von (14) 
endlich sei. diskutieren wir das Auftreten von Nullstellen des 
Reflexionsvermögens durch Aufsuchen der Nullstellen des Zählers 
g, sin?s + g, cos*s +g, sin2s=0. Der Zähler ist Null, wenn 

gleichzeitig sins = 0 und g,=0 oder gleichzeitig coss=0 und 
9 = (0 sind. Um evtl. außerdem vorhandene, nicht an zwei gleich- 
 geitige Voraussetzungen gebundene Nullstellen aufzusuchen, lösen 
- diese Gleichung nach sin?s auf und erhalten 


NER 


2 9." — 9%(9ı — 92) + 295 V9s" — 91 9 


sin? s = 
(9: — 92)" + 495" 


eer a Ein reelles s erfordert, abgesehen davon, daB der Ausdruck auf 
der rechten Seite zwischen 0 und 1 liegen muß, daß g,? > 9,9, ist. 
_ Wie aus der Betrachtung der ersten Näherung im folgendem hervor- 
gehen wird, tritt g, als ein gegen g,, 9, kleines Glied in Erscheinung, 
das im Rahmen der Gültigkeit unserer Näherung klein bleiben muß. 
Innerhalb ihres Gültigkeitsbereichs haben wir also keine weiteren 
Nullstellen des Reflexionsvermögens zu erwarten. 


§3. Vergleich der Näherungslösung mit der exakten Lösung 
für einen bekannten Spezialfall. 


M. Schlick!) und G. Bauer?) 


- A 
(17) n=—, 


für welchen die exakten Lösungen von (3) bekannt sind, das Reflexions- 
vermögen berechnet. Die Lösungen (5) sind in diesem Falle 


and 39. 1941 
1 
1 
(18) 
I 
I 
d 
p 
u 
t 
E 
lauf des Brechungsindex in III, nämlich 
— 
Py 
[223 
we 
1 V2 1 ) 
— + — — 
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wobei A eine Funktion der Schichtdicke d und der Brechungsindizes 
n, =n(z,) und n, = n(z,) an den Grenzen der inhomogenen Schicht ist 


Ng 
Ny — Ny 


Das Auftreten der Ww ug im Exponenten in (18) führt m einer 


Fallunterscheidung. Falls - ist, wird <A > und 


Ny Ne 
die Wurzel ist reell, so daß die nach (5) eingeführten Größen 
c @ 49 ce 1 
j= — |/ :— =—5-— Ww ver- 
und q(z) =. werden. Sie ve 


wandeln sich unter Wiedereinführen der WER (17)in p = < m = 


und q = z mit m= v5 . Die Größen (16) werden dann 


ho 
2 An A 


1 


Einsetzen der Größen (20) in (14') ergibt das exakte Reflexions- 


vermögen für den speziellen Verlauf (17) des Brechungsindex in der 


inhomogenen Schicht. 
. ale: In-n 
Diese Formel bleibt auch gültig, wenn 2 i ah. ist, 
wenn also 4 < + wird, Allerdings wird dann die Wurzel rein 
1 


imaginär und man führt statt m besser u = y- 


@ 

4 e 


benutzt (14°) in anderer Form 


21) + 9, TG? + 29,293 
Is + 9, s + 29, 


wobei 9,, Js, 93, 5 aus den Ausdrücken (20) hervorgehen, wenn man 
dort für m u setzt. Der Übergang vom ersten zum zweiten Fall 


| 
| 
19 
’ 
4 
6 
+ 
les € 
nn 
g, = (n,n, — n. ) 
en 2a A R au 
=(nn..+n one 
9, ( (NN, — Nyy 
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geschieht, wenn m = 3 =1 wird. Er tritt in 
Erscheinung, wenn man unter Festhaltung der Differenz n, —n 
der Brechungsindizes an den Grenzen des inhomogenen Mediums zu 
immer kleineren Schichtdicken iibergeht. Dann verschwindet der 
das Reflexionsvermögen modulierende Eintluß der Interferenz der an 
der Vorder- und Rückseite der inhomogenen Schicht reflektierten 
Wellen; gleichzeitig steigt das Reflexionsvermögen monoton an in- 
folge des Anwachsens von n’/n. 

Benutzt man statt der exakten Lösung in diesem Spezialfall 


die erste Näherung, so tritt an Stelle von p = — + -n der Nähe- 


rungswert p =, während q seinen Wert behält. — + > | ist 
aber gleichbedeutend mit rer —» 0, d. h. die Näherungsformel 


stimmt um so besser mit der exakten überein, je größer 7 wird, 


“0 


also bei fester Wellenlänge für hinreichend große Schichtdicken. 
> 1 geht, d.h. für 


ho 
4nA 


Die Näherung versagt, wenn 
i—e De 25, 
Ny Ng 
Diese kritische Schichtdicke D ist in normalen Fällen etwa eine 
Zehnerpotenz kleiner als 2,. Bei d = D geht die exakte Lösung in 
die Form (21) über. 


$4. Das Reflexionsvermögen bei einem beliebigen Verlauf 
des Brechungsindex in der inhomogenen Schicht 
Die in der exakten Lösung des Spezialfalls von $ 3 aufgetretene 
Fallunterscheidung ist von allgemeiner Bedeutung. u die Reihen- 
e 
o fy ~ 


Tr Tes . 
Unsere Näherung ist nur zu verwenden, wenn — fi <idh 


af 
=) <1, =) und wenn entsprechende Vorbehalte für 


entwicklung (6) verliert ihren Sinn, wenn 


dan n 
die höheren Näherungen, also die höheren Ableitungen von n ge- 
macht werden. Aber nach diesen Voraussetzungen gilt sie dann 
für beliebigen Verlauf von n in der inhomogenen Schicht. 
Wir betrachten im folgenden das Reflexionsvermögen in auf- 
einanderfolgend höheren Näherungen. 


A. Nullte Näherung 


In (6) sei nur f, von Bedeutung, alle höheren f;; (i = 1,2,3,...) 
seien zu In das (14) geh 
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über die Ausdrücke (16) nur die Werte von p und q ein, die sie 
an den Grenzen z, und z, der inhomogenen Schicht besitzen. Da- 
her kann der Brechungsindex im Innern der Schicht auch stärker 
veränderlich sein, wenn er nur an den Grenzen das in der ent- 
sprechenden Näherung verlangte Aussehen besitzt. In — 
Niherung ist p =m und q =0, d.h. in (16) ist zu setzen p, = n,, 


Ny = fras während q, = q, =0 d.h. n,’ 2 =0 anzu- 


nehmen ist. Dis Reflexionsvermögen (14’) wird dann 


var 


Ny N, — Nr N os- 8 + (nn — nn. 
(m 2 II iy (my 1 sin” 


22) 
(m N, + iy cos’ + (ny Nyy + My My)” sin? 8 


Man erkennt aus dieser Formel, daß in zwei Fällen das Reflexions- 
vermögen verschwindet, nämlich 


a) wenn sowohl — = "u als auch coss=0,d.h. 
n n 
I 2 


s= @)dz = 5 
> 


2 


oder der Lichtweg |n(z)dz= > (k +!/,) ist (vgl. Abb. 1). 


17 


Z | 
i Mi 


2 
Abb. 1. Ein Verlauf des Brechungs- Abb. 2. Ein zweiter Verlauf des Bre- 
index, bei dem das Reflexionsvermégen chungsindex, bei dem das Reflexions- 
verschwinden kann vermégen verschwinden kann 


b) wenn sowohl — = als auch sins=0d.h. 


oder der Lichtweg =k. * ist (vgl. Abb. 2). 
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Gl. (22) geht nebenbei für n, = ny = ny und s—>0 in die Fresnelsche 
Formel bei senkrechtem Einfall 


B. Erste Näherung 


Dabei berücksichtigen wir außer f, noch f, in der Entwicklung 
_ (6), so daß wir in die Ausdrücke (16) die Wertep=n und q= = a 

nn 
einzusetzen haben. Wir finden dann das Reflexionsvermögen für 
den Fall (vgl. Abb. 3), daß sowohl der Brechungsindex selbst als auch 


die ersten Ableitungen an den Grenzebenen springen. 


=! 


Abb. 3. Zwei Sprungstellen des Brechungsindex und seiner ersten Ableitung 


Die Koeffizienten (16) werden 


9 


= (MyM — + | 


= 94 — 4Ny Ny N2, 


(Ny Nz — + 
4n 


/ 


95 — Nı Ne, 


hy 


Wir erkennen aus diesen Ausdrücken, daß die Berücksichtigung der 
Sprünge der ersten Ableitung des Brechungsindex an den Grenz- 
ebenen nur Zusatzglieder der Größenordnung erster bis vierter 
Potenz von 2,/4a hervorrufen. In erster Potenz tritt 4,/4a nur 
in g, auf. Da das Auftreten von Nullstellen des Reflexionsvermögens 
nur von 9, bzw. g, abhängt, wird deren Lage also nur wenig von 
ihrer Lage bei der nullten Näherung verschieden sein. 


erh 
‘ 
[Ng 2.72 + 2n,n, 
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Aus den Ausdrücken (23) findet man für den Fall (vgl. Abb. 4), daß 
der Brechungsindex selbst stetig in die inhomogene Schicht übergeht 


n=", N2= Ny), dagegen die Sprünge in der ersten Ableitung 


erhalten bleiben (n,”+0, n,’ + 0) folgende vereinfachte Koeffizienten: 


= [ nin, + +(2 = 


2» 
(nin, — nim, ) =9g,— 
111 


Abb. 4. Zwei Sprungstellen der ersten Ableitung des Brechungsindex 


Wenn man in diesen Ausdrücken die jeweils in unserer Näherung 
überwiegenden Glieder berücksichtigt, wird das Reflexionsvermögen 


2 2 ..’\2 2 
(nin, + sin” 8 + (npn, — cos” 8 


+ nn, (nin, — ny n,) sin 28 


2 2 . . 
Der Faktor (2) in (24) bedeutet im allgemeinen eine starke af 
Verringerung des Reflexionsvermégens zweier Sprungstellen der tee ER 


ersten Ableitung des Brechungsindex gegenüber zwei Sprungstellen — 
des Brechungsindex selbst [vgl. Gl.(22)]). Das Reflexionsvermögen 2) 


verschwindet, wenn sowohl sins=0 als auch 


werden. 

Die Reflexion einer einfachen Sprungstelle der Ableitung erhält 
man in Übereinstimmung mit R. Gans?) aus (24), wenn man z.B. 
N,=0 wählt: 


(95) 


> 
9, 
= 9, 
dy 
ny 
r 
r A 
| 
"2 ” 4% 4 
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C. Höhere Näherungen 


Der Übergang zu höheren Näherungen ist durch die Entwick- 
lung (6) mit (7) zusammen vorgezeichnet. Z. B. wird in 2. Näherung 


26) (5) {3 (= 
ie und in 3. Näherung 


n 


[P@) = ne) + (2) 


4a 4a n n* 
ves .. 
In höheren Näherungen treten also fortlaufend höhere Potenzen von 
4,/4r in den Koeffizienten g, bis g, unserer Formel (14) für das 
Reflexionsvermögen auf. Außerdem treten von der 2. Näherung an 


Ef rea: Zusatzglieder gegeniiber s in (23) auf. 


7 Von Interesse dürfte noch das Reflexionsvermögen sein, das in 

7 2. Näherung durch zwei Sprungstellen des zweiten Difterential- 

 quotienten des Brechungsindex (n,”+ 0, n,”+ 0) allein hervorgerufen 

wird, während der Brechungsindex selbst und seine erste Ableitung 

stetig in die inhomogene Schicht einlaufen (nm, =, = My, 
n,= 0). Man erhält die Koeffizienten 


I ny Nyy 


und unter Vernachlässigung von Gliedern mit höheren Potenzen 
von für das Reflexionsvermögen 


+ sin” + (ning — cos 
4a 


4ny 

a i Dieses Reflexionsvermögen verschwindet, wenn sowohl sin s =0 als 


3 
n 
| werden. Das Reflexionsvermögen einer einfachen 
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Sprungstelle der zweiten Ableitung des Brechungsindex erhält man Be en 


daraus durch Nullsetzen von n. 


4 \4n 


Bei analoger Vernachlässigung von Gliedern mit höheren Potenzen 
von 4,/4% ergibt sich das Reflexionsvermögen einer inhomogenen 
Schicht, an deren Eingang und Ausgang sich je eine > 
des k.-Differentialquotienten des Brechungsindex 
\ dz* Jz=:, * 
befindet, sonst aber alle niedrigeren Differentialquotienten stetig 
übergehen, zu 


k n®))? sin? s + Cu ny - 8 


2° 


a 


4n 
und fiir eine einfache Sprungstelle des kten Differentialquotienten gilt 


4 \4n ny 


Fir ein inhomogenes Medium, dessen Brechungsindex selbst stetig 
und in allen Ableitungen stetig von n(—oo)=n; nach n(+ 00) = nq 
übergeht, gibt das angewandte Verfahren keine Reflexion. Diese tritt 
nur an Sprungstellen des Brechungsindex oder an Sprungstellen 
einer seiner Ableitungen auf. 


§5. Vergleich mit dem Reflexionsvermögen einer homogenen Platte, ; 
die in zwei homogene Medien anderer Brechungsindizes eingebettet ist 
Die allgemeine Reflexionsformel (14’) läßt sich auch folgender- 

massen schreiben 


4ny — Pi Pa = 
95 + (94 - 95) sin? s + 9g sin 28 


Eine homogene Schicht des Brechungsindex % zwischen zwei homo- 
genen Medien der Brechungsindizes n; und my hat das Reflexions- — 
vermögen 

4ny Nyy 


(34 =l-— 
) Q sin? 8 


entsteht aus dem allgemeineren Fall (33), wenn das eingelagerte 
Übergangsmedium zwar inhomogen ist, aber an den ite 
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nur eine Unstetigkeit der nullten Ableitung des Brechungsindex auf. 
weist; vgl. Abb. 1 und 2 


My, N, ~My, 


Dann ist 
4ny Nyy "2 


+ (nj ) n> 5) sin? 


Demnach stellt die von G. Bauer*) vorgenommene Vereinfachung 
seiner Gl. (12), die neuerdings von H. Schröder?) diskutiert wurde, 
den Übergang zu diesem Sonderfall dar. Der Übergang von n, 
nach n, in der inhomogenen Schicht muß außerdem hinreichend 
langsam erfolgen. Sonst wird das Glied mit sin 2s immer vor- 
handen sein. Es kann nur noch „zufällig“ verschwinden, wenn die 
Brechungsindizes und ihre Ableitungen an den Schichtgrenzen 
gerade so gewählt werden, daß g, in (16) Null wird. 


Herrn Prof. Dr. M. Czerny und Herrn Doz. Dr. H. Dänzer sind 
wir für ihre Ermunterung zu dieser Arbeit dankbar. 
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Die Berechnung des Feldverlaufs 
in einem ultrahochfrequenzdurchfluteten System 
konzentrischer Zylinder komplexer Leitfähigkeit 


Von R, Stachowiack 
eilung aus dem Kaiser-Wilhelm- Institut für Biophysik, Frankfurt a. M. 


(Mit 1 Abbildung) 
7 
Die Ermittlung der Feldverteilung in ultrahochfrequente Span 
nungen führenden Anordnungen ist für viele Fragen von Wichtig- 
keit. Sie ist gleichzeitig durch einen besonderen Umstand gegenüber 
ler gleichen Problemstellung in der Niederfrequenztechnik prinzipiell 
vickelter, nämlich dadurch, daß die mit dem Feld verbundene 
Strömung nicht rein dielektrischen oder Leitungscharakter hat, 
lern komplex ist, d. h. es entsteht eine Feldverteilung unter = 
chzeitiger Einwirkung des Leitungs- und des Verschiebungsstromes. i 
Neben physikalischen und technischen Fragen existieren ri 
allem in der Biophysik eine Reihe von Problemen, bei denen die _ 
Kenntnis solcher Feldverteilungen in ultrahochfrequenzdurchströmten 
Objekten komplexer Leitfähigkeit von Wichtigkeit ist. Es ist das 
der Problemkreis, der sich mit der Aufklärung des Wirkungs- 
mechanismus der Ultrakurzwellen auf den lebenden Organimus (Ultra - 
kurzwellentherapie) befaßt. Und zwar sind es sowohl makroskopische = 
wie auch Anordnungen in der Mikrostruktur des Gewebes, deren 
nähere Untersuchung hier von Interesse ist. Erstere sind gegeben 
in der Vielfachschichtung, die ein durchfluteter Körperabschnitt in. | 
seinem Aufbau aus Haut, Unterhautfett- und Muskelgewebe, Blut, 
Knochen usw. darstellt; letztere sind gegeben in der Zell- und Faser- — 
struktur des Gewebes. Für beide Fälle erscheint angesichts des 
verwickelten, geometrisch unregelmäßigen Aufbaues des biologischen Pe: 
Objektes die rechnerische Ermittlung solcher Feldverteilungen oe is 
allgemeinen sehr erschwert. Doch lassen sich immerhin — a 
einfache Anordnungen in grober Näherung auf geometrisch einfache 
Elemente zurückführen, die dann einer rechnerischen Behandlung 


zugänglich sind. Im folgenden ist eine solche rechnerische 


lung einer einfachen Anordnung durchgeführt, und zwar handelt es = 
sich um zwei lange konzentrische Zylinder im ursprünglich homogenen 


Feld, wobei drei verschiedene komplexe Leitfähigkeiten vorliegen. 
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7 Die Betrachtung knüpft an den Spezialfall des langen Voll. 
Zylinders an, den Krasny-Ergen') für Ultrahochfrequenz auf 
potential-theoretischem Wege ausgewertet und hierbei quantitativ 

die Möglichkeit des Auftretens von zusätzlichen Drehfeldern er. 
wiesen hat. 
Das innere Potential hat die Form 


>—ı 


bedingungen berechenbaren Konstanten F und @ enthalten die äußere 


be 1 normiert ist; r ist die Polarkoordinate. Die aus den Rand- 
und innere spezifische komplexe Leitfähigkeit k, bzw. k,: 


y ka — kj 

Die Komponenten der äußeren Feldstärke sind also im all- 
gemeinen nicht in gleicher Phase (zusätzliche Drehfelder). 

Im folgenden soll der Fall zweier konzentrischer langer Zylinder- 
körper im ursprünglich homogenen Feld, die theoretische Ersatz- 
möglichkeit durch einen homogenen Vollzylinder und zuletzt ein 
therapeutisches Anwendungsbeispiel gebracht werden. 


Wir haben hier zu unterscheiden (vgl. Abb. 1) zwischen einem 
Innenraum i, einem Zwischenraum z und einem Außenraum a; ¢, 


k, x und r sollen die oben angegebenen Bedeutungen haben. 
 £=r-cosa. 


Das Querschnittsverhältnis 9 = z ‚=1. 
Nun können in Analogie zum Vollzylinderproblem plausible 
Ansätze betrefis der Form der Potentiale gemacht werden. Die 


1) W. Krasny- ae HF-Techn. u. El.-Ak. 49. S. 195. 1937. 


| 
| 
jj 
dagegen ist das äußere | 3 
re Hierin ist x die Koordinate vom Mittelpunkt des Zylinders vom 
Radius R in Richtung der ursprüngli h 
d 
b 
N 
Pp 
d 
> 
Abb. 1. Schen nlezer Leitfähigkeit 
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Konstanten werden aus den Randbedingungen berechnet und zur 
Kontrolle überbestimmt. Die Ansätze sind folgende: 
Durch die konzentrische Einfügung des Innenzylinders in den 
äußeren Zylinder, der sonst ein homogenes Feld aufweisen würde, 
entsteht im Innenzylinder auch eines dieser Art: 
Dagegen kann das Zwischenfeld nicht allein aus einem neuen 
homogenen Feld bestehen, sondern hat ähnlich wie das Außenfeld 
eines einfachen Zylinders die Form er 
—9,=B-2+ 
Das äußere Feld ist formal sogar genau das. gleiche wie beim 


homogenen Vollzylinder 


Die Konstanten A, B, C und D dieser 3 Ansätze werden an 
den beiden Grenzradien R, und R, nun aus den Anschlußbedingungen 
bestimmt bzw. überbestimmt erstens für die Potentiale, zweitens für 
die senkrechten Stromkomponenten und drittens zur Kontrolle für die 
Parallelkomponenten der Feldstärken. Die ersten und dritten Be- 
dingungspaare sollen gleiche Ergebnisse bringen. 7 

P(R)=9.(R); 
. O9(R) _ 7 6 . 
Og (R) _ , (R) _ 
da 0a 

Durch Einführung der 3 y-Ansitze gewinnt man in der Tat 
fir A, B, C und D 4 Bestimmungsgleichungen und 2 Über- 
bestimmungen. Die erste und die dritte Zeile ergeben dabei, wie 
gefordert, die gleichen Bedingungen. 

A=B+C 
= k,(B—C) 
k,(B — C 9) = k, (1 — D) 

Durch einfache Auflösung erhält man | 

A = 4k k,/Nenner, 
B=2k,(k,+ k)/Nenner, 
= 2k,(k,— k)/Nenner, 

= (k, — k)-(k,+k)+%-(k, + k)-(k, kyNenner, 
+ k)-(k,+ k) + #(k, —k,)- 
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Aus den Formeln ist ersichtlich das elektrische Zusammen- 
wirken von 3 komplexen Leitfähigkeiten, die auf 2 konzentrische 
[2 


Zylinder und den Außenraum verteilt sind. 


ee In bezug auf die Wirkung im Außenraum kann der Doppel- 
zylinder durch einen homogenen Vollzylinder ersetzt werden, dessen 


äquivalente Leitfähigkeit k, sich durch Vergleich der 2 Konstanten D 
bezogen auf den allgemeinen und auf den speziellen Fall 


ka + ky 


leicht ermitteln läßt: 
k, =k 
711+ (1 Hk,’ 
ka hängt also nur von &, k, und k, ab. 

Als Anwendungsbeispiel aus der Ultrahochfrequenztherapie 
wählen wir den Röhrenknochen. Es interessiert uns Feld und Er- 
wärmung im Knochenmark und in der den Knochen weithin um- 
gebenden Muskulatur bei „ursprünglich“ homogenen Feld. Auch 
soll die äquivalente komplexe Leitfähigkeit des Knochens als ganzes 
in irgendeinem Außenraum betrachtet werden). 

Den Berechnungen sind die in der Tab. 1 aufgeführten elek- 
trischen Konstanten der Gewebe zugrunde gelegt. Setzt man diese 
Werte in die oben hergeleiteten Gleichungen ein, so ergeben sich 
die komplexen Feldkonstanten (Betrag und Phase). Für Muskulatur 
und Knochenmark sind diese, sowie die Nutzleistungen in Tab. 2 
aufgeführt. Es ist dabei für die ursprüngliche Feldstärke der Be- 
trag 1 V/cm angesetzt. 

Es ergibt sich also, daß, obwohl die Leitfähigkeit des Knochen- 
marks nur '/,, von der der Muskulatur ausmacht, die Erwärmung 
im Knochenmark nur auf etwa '/,, der der Muskulatur absinkt, ein 
Umstand, der seine Erklärung in der höheren, noch dazu nach 
quadratischem Gesetz wirksamen Feldstärke findet. 

Berechnet man auf Grund der hergeleiteten Formeln für einen 
homogenen Vollzylinder niederer DK. ohne Leitfähigkeit die Feld- 
konstante D, so ergibt sich ein mit dem in der Tab. 2 nahezu 
gleicher Wert (etwa 1). Damit ist bewiesen, daß die Wirkung des 
Röhrenknochens auf die Feldverteilung in der Muskulatur nahezu 
die gleiche ist, wie die eines isolierenden Vollknochens oder sonstigen 
sehr schlechten Leiters niederer DK. 


1) Das Problem ist idealisiert. Die natürliche Einlagerung des Réhren- 
knochens in die Muskulatur und insbesondere die üblichen Elektroden- 
anwendungen bringen Feldabweichungen. Weit distanzierte Schliephake- 
elektroden dürften den Idealfall einigermaßen verwirklichen. 
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Tabelle 1 
(f=50MHz und & = 0,5) 


Biologischer Muskulatur Knochenmantel | Knochenmark 


Körper 


Leitfähigkeit | 8,0 *) klein gegen 0,2 


in m S-em— 
DK. 90 **) etwa §8***) 


Leitfähigkeit 
2,5 etwa 0,22 


*, B. Rajewsky, Ergebnisse der biophysik. Forschung I, Thieme 
ipzig (1938). 

*) K. Osswald, HF-Techn. u. EI.-Ak. 49. S. 40. 1937. 

**) Landolt-Börnstein, Springer, Berlin S. 1043. 1923, Unter Ca-Ver- 
dungen. 


Tabelle 2 


Feldkonstanten | D (Muskulatur) A (Knochenmark) 


abs. Betrag 
Phasenwinkel 


Nutzleistung mW 


Als ein besonders eigenartiges Auswertungsergebnis findet man 
der Berechnung der äquivalenten Leitfähigkeit des Röhren- 
ochens 
kg = — 0,10 + 0,257 = + 0,27etir; (Key = 1139) 


in bezug auf einen Außenraum. Hierin ist die negative Verlust- 
komponente natürlich nur als rechnerische Größe zu betrachten, die 
sich nur auf das (äußere) Feld verzerrend auswirkt. Die Erwärmungen 
sind aber bekanntlich abhängig von den örtlichen as, Ohm- 


Zusammenfassung gy 


Zusammenfassend ist zu sagen, daB die obige Ableitung und 
das numerisch ausgewertete Beispiel zeigen, daß die theoretische 
Weiterführung des von Krassny-Ergen untersuchten Falles von 
in ein Medium eingebetteten homogenen Zylindern auf geschichtete 
Zylinder analytisch durchführbar ist und Aufschluß über die Energie- 
verteilung für ein praktisch wichtiges System zu liefern vermag. 
Es ergibt sich dabei, daß die Feldverteilung im Außenraum (Mus- 
kulatur) bei einem geschichteten Zylinder (Röhrenknochen) etwa die 
gleiche ist wie bei einem Vollzylinder. Die Nutzleist-ng im Knochen- 
mark ist nn der erhöhten inneren ne relativ groß. Die 
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ae oben behandelten ähnlichen Verhältnisse sind in mannigfacher 
Form überall bei der Ultrahochfrequenzdurchflutung biologischer 
Körper im Spiele, und zwar nicht nur in den makroskopischen Schich- 
_ tungen verschiedener Gewebe, sondern ebenso auch in der Mikrostruktur, 
z.B. in den schlecht leitenden Formelementen (Membranen) der einzelnen 
Fr Zelle. All diese Fälle sind allerdings aus den eingangs genannten 
_ Griinden einer rechnerischen Erfassung sehr schwer zugänglich. Hier 
besteht nur die Möglichkeit der direkten Messung, z. B. der unmittel- 
ait baren Erwärmungsmessung [vgl. etwa die Messungen von A. Stein- 
_ hiiuser’)], oder mittels Sonde in einem hinsichtlich geometrischer Ge- 
stalt, DK. und Leitfähigkeit identisch nachgebildeten Modell. Über 
ae solche Messungen ist a.a.O. [H.Schaefer und R. Stachowiak*) 
berichtet worden. 


1) A. Steinhäuser, Strahlentherapie 68. S. 79. 1940. 
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